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Abstract 

An overview of the momentum and frequency dependence of effective electron-electron inter- 
actions which favor electronic instability to a superconducting state in the angular-momentum 
channel I and the properties of the interactions which determine the magnitude of the temperature 
T c of the instability is provided. Both interactions induced through exchange of phonons as well 
as purely electronic fluctuations of spin density, charge density or current density are considered. 
Special attention is paid to the role of quantum critical fluctuations including pairing due to their 
virtual exchange as well as de-pairing due to inelastic scattering. 

In light of the above, empirical data and theory specific to phonon induced superconductivity, su- 
perfluidity in liquid He 3 , superconductivity in some of the heavy fermion compounds, in Cuprates, 
in pncitides and the valence skipping compound, is reviewed. To provide an anchor to the limits 
for T c , the solvable case of dilute fermions with interactions varying from very weak (BCS limit) 
to the unitarity scattering limit and beyond (Bose-Einstein condensation of molecules) realized in 
experiments on cold atoms in optical traps is also discussed. 

The physical basis for the following observation is provided: The universal ratio of s-wave T c 
to Fermi-energy for fermions at the unitarity limit for this last case is about 0.15, the ratio of 
the maximum T c to the typical phonon frequency in phonon induced s-wave superconductivity is 
of the same order; the ratio of p-wave T c to the renormalized Fermi-energy in liquid He 3 , a very 
strongly correlated Fermi-liquid near its melting pressure, is only O(10 -3 ); in the Cuprates and 
the heavy-fermions where d-wave superconductivity occurs in a region governed by a special class 
of quantum-critical fluctuations, this ratio rises to O(10~ 2 ). 

These discussions also suggest factors important for obtaining higher T c . Experiments and 
theoretical investigations are suggested to clarify the many unresolved issues. 

PACS numbers: 
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I. INTRODUCTION 

A. Scope of this Overview 

There are two distinct aspects to the theory of superconductivity [1] based on the Cooper- 
pair instability of the normal state of metals. The first is the theory for the kinetic energy and 
for the interaction vertices of electrons as a function of their momenta, spin and energy. The 
second is the solution of the model specified by the first aspect to understand and to calculate 
various properties of the superconducting state. In this overview, only the two simplest 
properties of the superconducting state are considered: the symmetry of the superconducting 
state and the transition temperature T c . They are the simplest properties since they are 
(in most cases) provided by the linearized version of the theory. I will discuss both the well 
understood problem of pairing vertices in the s-wave symmetry and transition temperatures 
due to electron-phonon interactions (e-ph) as well as the difficult and continuing problem of 
pairing due to electron-electron (e-e) interactions themselves in lower symmetries. 

The point of view taken here is that for both e-ph and e-e interactions, the second aspect, 
for every realized superconductor is solved by the Eliashberg extension [2] of the BCS theory 
to include the retarded nature of the effective interactions. The validity of this theory rests 
on the smallness of the parameter Xu c /W, where A is a dimensionless coupling constant, u c 
is the characteristic high frequency cut-off of the interactions and W is the unrenormalized 
electronic bandwidth. This parameter is well known to be small enough for the e-ph problem. 
I will argue that for every known case of superconductivity through e-e interactions as well, 
this parameter is small enough for the validity of the theory if not for accurate quantitative 
calculations. 

This leaves the first part of the problem, which is essentially the problem about the nor- 
mal state of the metal. The spin and momentum dependence of the interaction vertices 
can be projected onto the irreducible representations of the lattice. These together with 
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the single-particle spectral function A(k, uj) specify the possible symmetries of the super- 
conducting state. The strength of the vertices in the different irreducible representations 
and the frequency dependence of the vertices are required to address quantitative questions 
about T c and its relationship to normal state properties. Calculation of these quantities is a 
forbidding prospect. Yet systematic understanding of factors determining T c appears possi- 
ble with a judicious combination of empirical information, phenomenology and microscopic 
theory. 

One may adopt contrasting attitudes towards the question of T c . One attitude is that 
it is a foolhardy task born of naiivete. After all, we cannot even calculate the correlation 
energy of the normal state of a metal or of liquid He 3 , except by complicated numerical 
methods which are not our purview. Calculating the frequency and momentum dependence 
of vertices is even harder. The other attitude is that interesting issues arise in asking the 
question, and relating the answers to other measurable properties, and in examining the 
limit of validity of the answers. The second attitude is adopted here. This attitude is 
fostered by the enormous experimental and theoretical work on the e-ph problem on a wide 
variety of metals and the inter-relationship deduced in the various parameters determining 
their T c . In fact physical properties of the normal state of metals, for example the electronic 
effects on the phonon spectra, which depend on much more detailed understanding of e- 
ph interactions than T c , can be understood and calculated using controlled methods in 
agreement with experiments. The problem of e-e vertices is much harder. I will argue 
however that there is much to be learnt from the problem of e-ph interactions even for 
superconductivity due to e-e interactions. Also, there is a lot understood of the properties 
of the normal state of liquid He 3 through both measurements and calculations. We do have 
reliable deduction of a few Landau parameters for the whole range of Liquid He 3 densities and 
systematic attempts to calculate interaction vertices by constraining them by the measured 
Landau parameters and relating them to transport properties and the variations of T c with 
density [3] . For most metal of interest many more normal state properties are measurable in 
principle, and parameters fixing the essentials of the correlation functions of charge or spin 
or current correlations extracted. The quantitative question we shall be concerned about 
is the formulation of the model, the first aspect stated above, in terms of such quantities. 
Finally, it appears to be a fact that superconductivity of unconventional variety due to e-e 
interactions is present most often only in the region of the phase diagram of metals near 
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which a quantum-critical point [3] to some other phase lies, and a good empirical case can 
be made that it is promoted by quantum critical fluctuations. Critical regions, once their 
basic physics is understood, are always easier to understand quantitatively than regions away 
from criticality. This is because scale-invariance of critical fluctuations guarantees that far 
fewer parameters are required to specify the fluctuation spectra than away from criticality 
and these parameters can be deduced from experiments with less ambiguity than away from 
critical points. 

On the other hand, it is important to point out that calculations of the frequency and 
momentum dependent fluctuations for interacting fermions are fraught with difficulties. Es- 
pecially misleading are calculations of such quantities in perturbative approximations such 
as the random phase or self-consistent random phase approximations or phenomenological 
calculations which do not respect sum-rules. The inadequacies of such calculations are ev- 
ident for the case of superfluid He 3 and in cases where the calculations of simple models 
such as the Hubbard model with such methods are contrasted with those from elaborate 
numerical methods. The general moral is that analytical methods must be always tempered 
in this problem by the constraints of obtaining normal state properties in qualitative and 
quantitative form from the same fluctuations. 

There is a class of problems where rather precise numerical results on a model can be 
compared with the experiments and where at the unitarity limit of scattering universal 
results, independent of the details of the model are expected. This is the problem of low 
density of fermions on a lattice (or without a lattice) with variable inter-particle short-range 
attractive interaction in the s-wave channel. Definitive results obtained [5] on this problem 
bear comparison with experiments on cold atoms realized in optical traps [6]. The results 
serve as useful anchor to the discussion of some of the more difficult problems discussed 
here. 

After this prologue, let us look at the field of action. In fig. (JT]), measured T c for classes 
of superconductors is plotted on one axis and the measured specific heat coefficient 7 in 
the other. The black points refer to superconductors where T c is known primarily to be 
determined by e-ph interactions and the red points refer to the others. No direct correlation 
between the two quantities plotted is on the whole implied or discernible. The point of 
plotting 7 together with T c in the figure is simply a way of providing a parameter which is 
undoubtedly of importance for T c within a given class of materials as a judicious reading of 
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FIG. 1: Plot of Tc and specific heat coefficient for several classes of superconductors of interest. 
Black points denote those where the superconductivity is known to be promoted by electron-phonon 
interactions, red points denote those where it is likely that it is promoted by electron-electron 
interactions without the intermediary of phonons. 

this plot will reveal. In reading this figure, it should also be borne in mind that a prejudice 
has been made to include in the figure either those which the highest T c in their material class 
or those which will aid in the discussions in this paper. Included here are different classes of 
metals and metallic compounds including the historically first superconductor (Hg) and the 
old champion, the A15 compound (NbsSn, Ge) and the new champions, (Rb 3 Ceo and MgB 2 ) 
among those for which superconductivity is understood to be due to e-ph interactions. 
Among those shown in red include some cuprate superconductors, some heavy fermion 
superconductors, an Fe-pnictide as also the interesting case of doped BaBi0 3 . The crystal 
with the highest known T c m 164 kelvin, H gBa2Ca m -iCu m 02m+2+i (Hg l:2:m-l:m) with 
m=l, 2, and 3, under pressure of about 50 GPa |7], is not in the plot; its specfic heat at 
such pressures appears not to have been measured. Superffuid liquid He 3 , not in the plot 
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but which will be discussed, has a T c pa 2AmK near the melting pressure with 7 of about 
50Joule/moleK 2 . Dividing the latter by the ratio of the He 3 mass to the electronic mass, 
the equivalent electronic 7, is about 6mJ/moleK 2 . 

The observed maximum T c from e-ph interactions, realized in metals like Pb — Bi alloys, 
NbsSn, MgB 2 , etc., is of order 6^/10. Here 6d is the Debye temperature. For the e-e 
promoted cases, T c , for both the heavy fermions and the Cuprates is 0(E / /10 2 ) , where Ef 
is the renormalized Fermi-energy. For He 3 , the same ratio is about Ef/10 3 . We will discuss 
that there are good reasons for these orders of magnitudes. I will also argue that He 3 is 
about as high a value of T c /Ef as one is likely to get from fermion interactions unless one 
is at proximity to a quantum critical point of a special kind, which is the case for both the 
Cuprates and the heavy-fermion compounds. 

B. Organization and Summary of this overview 

I will first present in Sec. II, the minimal necessary summary of the theoretical back- 
ground of aspects relevant to pairing symmetry and T c . I will start with the general vertex 
in the pairing channel and discuss the basis for its approximation by the form generally 
used. I will then present the parameters that go into determining T c through Eliashberg 
equations for general forms of effective electron-electron interaction vertices. The param- 
eters determining T c for the electron-phonon problem and in the general case will also be 
summarized. The favored "angular-momentum" state I (more properly the irreducible rep- 
resentation of the point group) of pairing will be extracted from the momentum and the 
spin-dependence of the interaction vertices. Also discussed here, following Kohn and Lut- 
tinger [S], is the instructive case of superconductivity parameters in jellium for pairing in 
any angular momentum. 

This, the determination of the favored channel of pairing is the easy part of the prob- 
lem. The hard part is the discussion of the factors determining T c . We will review that 
they are the spectral strength of the fluctuations exchanged, their characteristic energy and 
their distribution with respect to thermal energies, the coupling of these fluctuations to 
fermions and the density of states of the latter over the scale of the characteristic energy of 
the fluctuations. One of the important points emphasized here is that these factors are not 
independent. This is explicitly evident in the empirical evidence presented (and its deriva- 
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tions given in the Appendix) for the e-ph induced superconductivity. The same is expected 
for e-e induced superconductivity. 

I will discuss the variation of T c with the variation in the distribution in frequency of 
the spectral weight of the fluctuations, subject to their integral being held a constant. This 
includes the discussion of the effects of inelastic scattering in depressing T c (and raising the 
ratio (A/T c )). Important differences between s-wave and finite angular momentum pairing 
is evident from this analysis. Using this it will be shown that the simplest class of quantum- 
critical fluctuations (called Gaussian fluctuations here) are bad for T c . In this class the 
frequency spectrum of the fluctuations is peaked near thermal frequencies. I will review 
that the form of critical fluctuation spectrum consistent with the normal state experiments 
in cuprates and the heavy-fermions and quite likely the pnictides does not suffer from this 
limitation. 

Armed with this general knowledge, I pass on to the semi-empirical information on vari- 
ous classes of materials and its understanding based on the above. First superconductivity in 
transition metals and compounds through e-ph interaction is reviewed, in Sec. III. There are 
several things to be learnt here because a variety of experiments provide inter-relationships 
between parameters which determined T c . Simple theoretical ideas have been used to un- 
derstand such relationships. I argue that this understanding is useful for the e-e problems 
as well. 

In preparation for pairing in metals due to e-e interactions, I will review the connection 
between vertices of Fermi-liquid theory, in particular the Landau parameters and the pairing 
vertex in Sec. IV. This will be accompanied by empirical information and calculations on 
the superfluidity of He 3 . I will then discuss superconductivity through other interactions. 
The discussion will be based both on data, empirical considerations and the theoretical 
guidance which is available. Superconductivity in Heavy-Fermions, Cuprates, Pnictides and 
the interesting case of Bai_ x K x Bi0 3 , which appears to have electronically induced s— wave 
pairing, are successively discussed in light of the theoretical discussion. Major unsolved 
experimental and theoretical issues are highlighted. 

What limits the T c /Ef for finite i e-e induced pairing? There is first of all the deleterious 
effect on T c of the normal self-energy which is principally determined by the coupling constant 
A in the s— wave channel. Generally this is larger than the coupling constant for the pairing 
self-energy A^. This is because the important effective interactions are always short-range. 



9 



Second the collective modes generally have only a fraction of the total spectra weight. Third, 
due to the part of the excitation spectra at thermal frequency, the role of inelastic scattering 
due to real as opposed to virtual scattering in pair-breaking is stronger for finite i pairing. 
The gain in T c from e-e induced superconductivity of a larger characteristic energy giving a 
large prefactor is thus mitigated and special conditions are required for a substantive T c . 

The empirical results on e-ph promoted superconductors, and on classes of e-e super- 
conductors: liquid He 3 , Heavy-fermions, Cuprates, Pnictides and Valence- Skippers is re- 
viewed in light of the discussions summarized above. I will explain how we understand 
semi-quantitatively that the maximum T c is only an order of magnitude lower than the 
cut-off frequency for e-ph superconductors. This is also about the upper limit obtained for 
the attractive Hubbard model with the same cut-off. For liquid He 3 , the lower value T c /Ef 
of O(10~ 3 is understood from the coupling constants obtained from the appropriate gener- 
alization of measured Landau parameters as due to the larger reduction from self-energy, 
the reduction of the cut-off frequency by the large renormalizations and the fact that the 
collective fluctuations have only a fraction of the total weight of the fluctuations given by 
the sum-rules. 

We gain a factor of O(10) in this dimensionless ratio both for the heavy-fermion super- 
conductors as well as the Cuprates and nearly that for the Pnictides. This is undoubtedly 
related to the empirical results emphasized in this overview that that high T c in Cuprates, 
in Heavy- Fermions, the Pnictides and the Valence-Skipper is related to proximity to a qcp 
of unusual variety. As noted, quantum-criticality of the Gaussian kind is deleterious for 
T c because it reduces the cut-off in the fluctuations and promotes pair-breaking through 
inelastic scattering. The unusual nature of the fluctuations responsible for the relatively 
high T c of Cuprates and (in dimensionless terms) of the heavy-fermions is that the spectra 
though critical is distributed over a wide frequency range and that it is nearly momentum 
independent. The former leads to lower inelastic scattering as temperature decreases while 
the latter gives the smallest ratio of I = to finite I coupling constants. For Cuprates direct 
evidence also shows that the cut-off frequency of this unusual spectra is only a factor of 
about 4 below the fermi-energy. 

The derivation of the spectra for the relevant qcp for the Cuprates [9j [10] relies on the 
discovery of an unusual competing order parameter in the Curpates. It is shown that the 
Action for the model can be written as a product of orthogonal topological variables, one 
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of which has local interactions in time and logarithmic interactions in space and the other 
which has local interactions in space and logarithmic interaction in time. The general class 
of microscopic models where such properties hold for the qcf is not known; it appears em- 
pirically to contains the microscopic models that describe the Heavy-fermions and possibly 
the Pnictides and the valence skippers. This is an exciting new development in the study of 
critical fluctuations which needs further thought and work, both experimental and theoret- 
ical. 

The case of the valence skipper Ba x Ki_ x BiO^ is discussed as it appears to be an electron- 
ically induced s-wave superconductor. Valence skippers have e-e induced pre-formed s-wave 
Cooper pairs in the normal state. Superconductivity consists in obtaining phase coherence 
of such pairs. There is some evidence, for which there is no theoretical understanding, that 
their fluctuation spectra is similar to that in the Cuprates. If so, not only is inelastic scat- 
tering not very deleterious, the normal self-energy from virtual processes is not as hurtful 
as in finite i pairing. To my mind discovery of other valence skippers with higher electronic 
densities hold the best promise of further increase of T c . 

The question is often asked, is room temperature superconductivity possible? The fact 
that we are only a factor of 2 away in Cuprates and we have not progressed beyond that 
in the past 15 years is both promising and depressing. My answer however is yes, and 
that the best prospects for reasons given in this paper are for electronically induced s-wave 
superconductors. The upper limit for T c /Ef is provided by the results on the attractive 
dilute Hubbard model near unitarity, which as reviewed here is ~ 0.15. I can only hope 
that, for purposes of large scale applications, this will happen in material classes which are 
three-dimensional, malleable in their bulk form and easy to fabricate. 

II. THEORY OF PAIRING SYMMETRIES AND OF T c 

Soon after BCS theory, Eliashberg [2] used the field theory methodology developed for 
superconductivity by Gorkov [11 J to formulate the theory of superconductivity to include 
the frequency dependence of the effective interactions through exchange of phonons. This 
theory is valid if (Xu c /W) << 1. The physical content of this limit is the Migdal theorem 
which proves that the vertex corrections in the theory of electron-phonon interactions are 
of 0(Xu c /W) compared to the leading dimensionless vertex A. 
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The experimental proof that the superconductivity in metals such as Pb, Sn, Nb, etc. is 
induced by electron-phonon interaction rests on analysis of tunneling spectrum [T21 H3] in 
these metals using the Eliashberg theory and the measurement of the spectrum of phonons 
by neutron scattering. The theory also provides experimental proofs of its limit of validity. 

With one important modification, which does not affect the linearized Eliashberg theory 
which is enough to determine T c , the theory can be used for pairing in any symmetry of 
degenerate fermions due to exchange of any kind of fluctuations, provided they can be 
regarded as bosons, and provided the Migdal condition is satisfied, where oj c is the upper 
cutoff of the spectrum of bosons. 

Excellent reviews of the technical aspects of the Eliashberg theory have been written [H] . 
To fulfill the limited goals of this review, I need to present only the dependence of some of 
the results of the theory on the parameters of the starting model. The model specifies a 
spectrum of fermions, a spectrum of bosons, and the interactions between the fermions and 
the bosons: 

H f = e k C ka C k- (1) 

k.cr 

The transition temperature is the temperature at which the particle-particle scattering 
vertex with total momentum and energy diverges. The exact vertex Ts{k, k + q) for 
fermions with total spin S = 0,1, scattering from k, — k to k + q, — k — q, with energy u 
to u) + v follows the Bethe-Salpeter equation through which it is related to the irreducible 
particle-particle vertex Is(k, k + q), as shown in the Fig. ([2]): 




FIG. 2: Bethe-Salpeter Equation for the relation between the complete particle-particle vertex V 
and the irreducible vertex /. 



T s (k, k + q) = l s (k, k + q)- (T/2) ]T I s (k, p)G(p, u)G(-p, -u)T s {p, k + q). (2) 

p,uj 
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Here Is ( k,p) is the irreducible vertex in the particle-particle channel with total momentum 
zero, which means it contains all diagrams which cannot be cut into two parts by cutting 
two particle lines (right-ward going lines in the convention of Fig. (|2])). G(p, lo) is the 
single-particle Green function at momentum-energy p,o;. In ([2]), and fig.pj), the integral 
over energy has a cut-off energy of the vertex Is- This cut-off u c is the upper energy scale 
of the collective fluctuations of the fermions or of other modes with which the fermions 
interact. For problems of interest, we expect this cut-off energy to be much smaller than Ef 
but much larger than T c . This can in general only be justified aposteriori. 

Quite generally, for the e-ph or the e-e problem, one can define a four-fermion inter- 
action vertex Is(k, k + q) which is irreducible in the particle-particle channel through the 
interaction Hamiltonian 

H int = 1/2 ^29aA k > k + qK,<5(- k ' ~ k - 0)^/3,7,5 (q, ^ n )c+ +q ^C+ k+q7 C_ M C kiQ 

(3) 

•^a,/3,7,(5(q) w n) is the propagator of the fluctuations which are exchanged by the fermions, co n 
sums over the boson Matsubara frequencies, and g at p{\t, k + q, v) is the interaction vertex 
with v denoting any additional idex, such as the polarization for phonons, needed to specify 
the coupling to the fluctuations. For interactions with phonons or any other density or 
current fluctuations a = /3;j = 5. For isotropic interactions with spin-fluctuations or 
spin-current fluctuations, the spin-dependence in the fermion operators and the Fluctuation 
propagator in has the form cr 7 < 5 ■ crg iCf . 

In general three integral equations specify the Eliashberg solution: the equation for the 
normal self-energy of fermions, the equation for the pairing self-energy of the fermions, 
and the equation for the self-energy of the Bosons due to superconductivity. The latter is 
unimportant for the e-ph problem because the corrections due to superconductivity are of 
order (m/M) 1 / 2 . But they are likely to be very important for e-e problem below T c , because 
superconductivity itself changes the spectrum of collective e-e fluctuations by gapping the 
low energy single-particle spectrum. This will have quite significant effects, for example on 
the temperature dependence of the superconducting gap. However, since we are interested 
here only in T c and that too only in a mean-field approximation, this will be neglected. 
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A. Considerations on Pairing Symmetry 

1. Pairing Symmetry due to e-ph Interactions 

Let us start with the familiar example of electron-phonon interactions when the irre- 
ducible interactions are specified by 

H e . ph = J2 ^/(ZMu^I^c+^cUKv + &± q „), (4) 

frqio bq U are the phonon creation and annihilation operators and (^(k, k + q) = 
yjh/ (2Mw qv )^ is the electron-phonon vertex, M is the ion-mass and w qi/ are the frequen- 
cies of phonons of momentum q and polarization v. A specific form for the electron-phonon 
vertex will be introduced later. 

Optical phonons are to a first approximation dispersion free. Their interactions with elec- 
trons in a metal are effectively 5-function in real space, which only gives s-wave scattering. 
The interaction with acoustic phonons must vanish in the long wavelength limit. Moreover, 
the partial density of states of phonons is peaked near the zone-boundary. The effective 
interaction is therefore short-range, of the order of the lattice constant. As will be discussed 
later, this is true for pseudo-potential metals in which the e-ph interactions are weak. For 
transition metals and compounds and other co-valent metals with larger interactions and 
larger T c , the e-ph interactions are even shorter-range. Dominant attraction is therefore 
always in the s— wave channel. 



2. Effective electron-electron repulsion in the s-wave channel in Jellium 

Electrons also repel. The effective dimensionless electron-electron interaction (normalized 
to the density of states at the chemical potential), and assumed to be cut-off in energy only 
at the Fermi-energy, 

H = N(0) < 47re 2 /(g 2 e(g)) > (5) 
where e(q) is the dielectric function and the average is over initial and final states at the 



chemical potential in the same manner as in the definition of the A in Eq. (24). It is 



interesting to consider the relationship of \x to A with both calculated for jellium. The 



dimensionless electron-phonon coupling constant A for jellium, defined in analogy with (23) 
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is 

\ e ~ ph = N(0) < g 2 scr {q) > /M < oo(q) 2 >, (6) 

where g SC r{<l) is the screened electron-phonon coupling function, related to the bare el-ph 
interaction function g(q) by 

9scr(q) = g(q)/e(q), (7) 
and the average is over the same states as in[5j For jellium, g(q) is given by 

g 2 {q)/{Mnl) = 4vre 2 /g 2 , (8) 

where the bare ion plasma frequency Q p is related to the actual frequency of the ions in 
jellium by the Bohm-Staver relation, 

u\q) = n 2 p /e(q). (9) 

We therefore have [15] 

\ e - ph = fx. (10) 

In view of the fact that jellium is completely characterized by one parameter r s , such a 
relation should not be surprising. One could go further and say that, this relation could be 
projected onto any angular momentum pairing channel I of electrons scattering from k, — k 
to k', — k' and therefore that such a relation exists in each channel, X e e ~ ph ~ fj,g. For Fermi- 
Thomas screening e(q) oc q 2 / (q 2 + q^p), where q^ F is the inverse screening length. One may 
project the interactions on to pairing channel I and find that the effective dimensionless 
direct repulsion as well as e-ph induced attraction parameters both are proportional to 
exp(-^). 

Cohen and Anderson [16], in a not easily accessible paper, which has several interesting 
ideas, have discussed how in a lattice, Umklapp scattering which in effect is due to the 
modulation of the charge density within a unit-cell due to phonons or local field effects, 
increases the coupling constant \ e ~ ph relative to fi. This paper emphasizes that strong e-ph 
interactions mean the modulation of electronic bonds. This point will be discussed further 
in Sec. (Ill) in the context of transition metals and compounds, where this idea is most 
prominently borne out in the empirical data and implemented efficiently through the tight- 
binding representation of e-ph interactions to provide a theory of e-ph interactions as well 
as a quantitative theory of the anomalies in the phonon spectra [17, 18J in transition metal 
and compounds. 
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We should also note the important fact that when the effective e-e interactions in the 
pairing channels are expressed so that they are retarded only over the same range, < 00 >, 
as the e-ph interactions they are reduced to 

// =fi/(l + fihi(Ef/ <u>). (11) 

3. Pairing in I 7^ through Incoherent particle-hole fluctuation exchange 

In an interesting paper Kohn and Luttinger [8] showed that the effective electron-electron 
interactions must, due to the sharpness of the fermi-surface, always be oscillatory in the 
momentum transfer 2/cpcos(#) where cos(#) = kf ■ k f , with both k and k' are close to the 
fermi-surface. From this it follows that there is always attraction in the pairing channel at 
some or the other angular momentum angular momentum £. The issue then is whether it is 
of sufficient magnitude to give T c larger than that from typical e-ph interactions. 

Kohn and Luttinger also provide an estimate for the effective interaction when the bare 
interaction is weak. The direct (screened for e-e interactions) repulsive interaction falls of 
exponentially in £, while the oscillatory exchange attractive interactions fall off only as l/£ 4 , 
the latter therefore wins for large enough £. The same is true for hard-core interaction 
of finite radius. For He 3 , taking this radius equal to the diameter of the He atoms, £ = 1 
channel is already attractive and provides using a BCS type expression, a T c »s o; c exp(— 2.5). 
Here u c is the cut-off. If one takes u c ~ Ef, one gets a factor 10 2 too large a value compared 
with the experiments. There are several things wrong with this. For £ 7^ pairing, one 
must include the repulsion in the £ = channel in the normal self-energy, even in the 
weak-coupling limit. This cuts down the estimate by about two orders of magnitude. The 
leading effect over the weak-coupling limit is that, the quasi-particle renormalization, due to 
increasing interactions, also renormalizes the attractive coupling constants in the exponent 
downwards as well as the cut-off cu c . 

An important point is that for any reasonable interactions, the Kohn-Luttinger result 
that the attractive interaction falls off as l/£ 4 is unlikely to change much in more com- 
plicated calculations with incoherent particle-hole interactions. Then, for example, for the 
parameters of He 3 , the £ = 2 pairing has a transition temperature of about 10~ l7 K, even 
without considering the s-wave repulsion. One may conclude that incoherent particle-hole 
fluctuations are not a very good way to get to high T c . How then to make use of the 
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higher cut-off energy of electronic excitations. The alternative to incoherent particle-hole 
fluctuations are collective modes of electronic fluctuations, to which we turn next. 

4- Exchange of collective fluctuations 

Part of the spectral weight of particle-hole fluctuations appears in collective modes if 
the effective interactions are strong enough compared to the kinetic energy. In general a 
particle hole-fluctuation is characterized by an internal momentum k and a center of mass 
momentum q, besides other indices which specify the channel of the excitations, for example, 
density or spin or current or spin-current or inter-valley, interband, etc. A useful limit is 
obtained from the fact that a bound state of a fermion particle and hole (and indeed particle 
and particle or hole and hole) has Boson commutation relations. In a bound state, there 
is no dispersion in energy at any q as a function of k, so the latter index may be summed 
over. To regard the fluctuations as collective rather than incoherent, one must make sure 
that the dependence on k is unimportant to do sensible calculations. Also important is the 
integrated spectra of the collective modes and the incoherent particle-hole spectral weight of 
fermions for any given channel, is generally fixed by sum-rules. It is sinful double- counting to 
introduce collective variables, for example spin-fluctuations, and give them the total weight 
of S(S + 1) per spin S, and do calculations of their interactions with fermions in the same 
way that one does for phonons. 

5. Symmetry of Pairing Interaction from Exchange of Spin-Fluctuations 

As we will see later, the frequency dependence of the interaction (with total spectral 
weight and coupling constant fixed) is of much more importance in determining T c in finite 
I than for i = 0. The anisotropy of the interaction as well as the details of the geometry of 
the fermi-surface(s) are also important in determining the irreducible representation favored 
for pairing. However, essential aspects of finite i pairing, specifically the distinction between 
the pairing promoted by ferromagnetic and anti-ferromagnetic fluctuations is revealed in a 
simple model calculation [19] based on a frequency independent and isotropic interaction in 
a model of a spherical fermi-surface. Consider effective electron-electron interactions of the 
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simple form: 

Him = 1/2 J(|k - k'DfT^ • 0- 7 5 C k+q/2,a C -k+q/2,7 C -k'+q/2,<5Ck'+q/2,/3 (12) 

kk'q 

In the weak-coupling approximation, the partial wave components of the pairing interaction 



Vg follow from Eq.(|12|) to be 

3 j 2^ dxxP £ (l -2x 2 )[-J(2k F x)}, (13) 

where the pre-factor is 3 for spin-singlet (even £) pairing and —1 for triplet-spin (odd £) 
pairing. 




One may analyze the consequences of the q- dependence of J(q) in of Eq.(13) using 
fig. (|3] ). The ferro-magnetic interactions have a peak of —J(q) at q = 0, while the anti- 
ferromagnetic interactions have a peak at the zone-boundary. First note that if J(2kpx) is 
independent of x, Ve = for all I. This is because J(q) independent of q implies a delta- 
function interaction in real space where any finite i Cooper pair wave-function has zero 
magnitude. On the other hand for i = 0, a constant J(q) represents a magnetic impurity, 



which suppresses pairing. One also can see from Eq.(13) and the smooth variation of J(q) 
in fig. ([3]) that Vo > 0, V\ < and V2 > for ferro-magnetic interactions. In other words 
both isotropic and anisotropic singlet states are disfavored by such interactions, while triplet 
i = 1 state is favored. For anti-ferromagnetic interactions, Vo > 0, V\ > 0; both isotropic 
singlet and the triplet state is disfavored. Whether V? < or > 0, depends on the detailed 
form of J(q). A strong enough peak in J(q) near the zone-boundary favors such pairing. 

It is straightforward to extend this analysis [TS] to the more realistic case for the momen- 
tum dependence of the interaction taking the crystal symmetry into account. It is expected 
that the antiferromagnetic fluctuations peak at a point q = Q along a symmetry axis and 
so does J(q). I refer to the original paper for some details. An important point to emerge 
from that analysis is that the pairing symmetry depends both on how steep is the increase 
in J(q) near the zone-boundary and on the details of the electronic dispersion e(k) near 
the chemical potential. For example the electronic structure may choose between extended 
- s or d-type pairing. This point appears relevant to the case of the recently discovered 
Fe-pnictide superconductors, where the band-structure with hole-pockets at the zone center 
and electron pockets near zone corners may favor nodeless pairing with a phase difference of 
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FIG. 3: q-dependece of FM and AFM interactions 



7r between the gaps at zone-centers and the gaps at zone corners. It appears possible that 
for appropriate band-structure AFM fluctuations may even promote triplet pairing. This 
may be implicated in the triplet superconductivity in SV2-RWO4 [20], which shows no FM 
correlations but modest AFM correlations. 

One can physically understand for the simple square lattice band-structure why a steep 
increase of J(q) near the zone-boundary favors d-wave. The wave-function for this interac- 
tion favors anti-parallel correlations of spins on nearest neighbors. Such a spin-correlation 
is also produced by the d— wave BCS wave-function, for such a band-sturcture. 

A generalization of magnetic fluctuations due to spin-moment correlations to magnetic 
fluctuations due to orbital-moment correlations will be discussed in connection with the 
physics of the cuprates in Sec. (IV). 

While the arguments above may provide valid grounds for discussing pairing symmetry, 
they are no help in thinking about the value of T c . For that one must turn, as we do next, 
to the frequency dependence of the pairing vertex, which in most cases of interest is given 
by the frequency dependence of the collective modes exchanged by the fermions. 
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B. Effect of the Frequency Dependence of Fluctuations on T c in s-wave and higher 
angular momentum pairing 



We know that we do need for high T c : 

(A) Large spectral weight in the collective fluctuations being exchanged by the fermions, 

(B) Large coupling of these fluctuations to the fermions, 

(C) Large frequency scale of the fluctuations and 

(D) Large density of states of fermions around the chemical potential in a range of 
frequencies of the upper cut-off of the fluctuations. 

(A) is fixed for phonons but not for e-e fluctuations and must be carefully considered so as 
not to overcount the degrees of freedom and for its effect in renormalizing (B). We will see in 
Sec. Ill that (B), (C) and (D) are not independent in the e-ph problem. Arguments will be 
given that this is true also for the e-e problem. For the case of e-e interactions there is also the 
additional and important consideration of the distribution of the spectrum of fluctuations, 
due to the role of inelastic scattering and the difference of normal and anomalous self-energy. 
We proceed to discuss this immediately. 

To discuss the role of the spectrum of fluctuations one must turn to the careful analysis 
of the solutions to the Eliashberg equations. This was done in a beautiful paper for s-wave 
pairing through phonons by Bergmann and Rainer [21], followed up to examine the role of 
AFM fluctuations in s-wave pairing [22] arid finally for d-wave pairing [23]. 

Let us start with the simplest situation when the spectrum of fluctuations is extended 
more or less uniformly over a large range up to a cut-off u c » T c . A simple generalization 
[23] 124] of the McMillan approximate solution [25] [2S] to £ ^ pairing gives 



for A S ,A^ of 0(1) or smaller. A^ is the dimensionless coupling constant in the £-th particle- 
particle channel in which pairing is presumed to occur and which therefore occurs in the 
pairing self-energy in the Eliashberg equation. Since the self-energy must respect the full 
symmetry of the lattice this always includes the I = channel. For a square lattice in two 
dimensions or a cubic lattice in three dimensions, the first £ ^ momentum channel included 
in A s is the fourth. Since in any reasonable model, the coupling constants go down fast with 




-1 + A 



(14) 
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we need only to include the I = in it. Eq.(14) has the consequence that for short-range 



interactions, T c /u c is smaller for £ ^ than for I = 0, with the decrease depending on the 
range of interactions. 



Just as Eq. ( 14 ) does not include the effect of pair-breaking due to any magnetic impurities 
present for i = wave or due to any form of impurity for I ^ 0, it does not include the effect 
of pair-breaking due to inelastic (real processes as opposed to virtual processes) scattering. 
These processes weigh the spectrum of fluctuations over the range that their frequency is 
smaller than 0(T). More detailed analysis of the Eliashberg equations is required to study 
this. 

A quantity of interest to study the role of the frequency dependence of the fluctuations 
is the functional derivative of T c with respect to the change in the spectral function A{uS) 
[21] which appears in the kernel of the Eliashberg integral equations: 

u5A(u) (15) 

For s-wave pairing, the analysis of Bergmann and Rainer shows that > for any 

frequency for s-wave pairing, reducing to in the limit w — > 0, consistent with Anderson's 
theorem about the effect of weak impurity scattering on T c . Increase of spectral weight in 
any region of the frequency spectrum increases T c for s-wave pairing. The maximum value 
of this quantity appears to occur near u ~ 2nT c . 

For finite £, the situation is different and depends on the ratio g = A^/A s . Then, for 
g < 1, one can define a frequency 



cu*^T c ex V (l/g) (16) 

such that for to > to*, Wjf^y > 0, while for to < uj*, w^f^y ^ 0. The physical reason 
for this is that thermal occupation of the fluctuations (other than those of long wavelenth 
which only lead to forward scattering) has similar effect as impurity scattering. This is 
pair-breaking for I ^ because it mixes the phase of gap at different parts of k-space in 
the pure limit. These conclusion are found in the solution of the Eliashberg equations for 
an Einstein model of spin-fluctuations, see Fig.Q. Note, also the much larger reduction in 
T c compared to the zero-temperature superconducting gap A(0), leading to a large value of 
A(0)/T c characteristic of the Cuprates. 
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FIG. 4: Dependence of T c /T c q and A/Ao on inelastic scattering for the case of d-wave pairing. 
An Einstein spectrum of spin-nuctutations at uje is used for simplicity with frequencies 1 and 
additional Einstein modes are introduced with frequencies (in units of we) as marked for curves 
a, b and c. T c o, and Ao are the values without inelastic scattering and so is the unmarked curve. 
Dashed and full lines represent approximate and exact solutions of the Eliashberg equations. For 
details, see ([23]) 

C. Quantum-Criticality in Relation to T c in e-e Induced Pairing 

The highest T c 's in the Cuprates are undoubtedly through electronic fluctuations, and 
near a quantum-critical point (qcp), as are those in the heavy-fermions. This as we will 
see is likely be true for the pncitides as well. As discussed below for each of these cases, 
we know this because the observed transport and thermodynamic properties for T ^ T c in 
these compounds can only be understood as due to scattering of fermions from fluctuations 
which have a singularity in the limit u> — > 0, T — > 0. It is therefore important to ask about 
the role of spectral distribution of fluctuations near qcp's in determining T c in light of what 
we have learnt in the previous section. 

Much less is understood about the universality classes of quantum critical points (qcp) 
than classical critical points where the scale-invariant frequency and momentum dependence 
of critical fluctuations for the classical critical points has been catalogued into different Uni- 
versality classes [27] . A qcp is the point p = p c , where by varying a parameter p, for example 
pressure or electron density, the transition temperature to some broken symmetry — > 0. A 
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simple generalization of the dynamics of classical critical phenomena to the quantum case 
[U [2HH3I] (which may be termed a Gaussian qcp since the fluctuation spectra is determined 
by the renormalized coefficient of the quadratic term of the variable describing the fluctuat- 
ing order parameter) is as follows: The spectral function of the fluctuations near a qcp may 
in general be be written in the scaling form: 

X n (q,u>,p,T) =XoC dAI J r i( k ZM z ,u/T), (17) 

£ r is the correlation length at T — > which diverges as 

Z r <x\ P - Pc \-», (18) 

and k is measured from the wave- vector of the symmetry breaking, du = d, the physical 
dimension if the order parameter is that of a conserved quantity, otherwise it can be differ- 
ent. The correlation length in the time-direction or equivalently the frequency scale of the 
fluctuations is given by 

& oc £. (19) 

z is the dynamical critical exponent which scales the frequency of fluctuations to their 
characteristic spatial extent reflecting that the dispersion of the critical fluctuations depends 



on the nature of the broken symmetry and in general has the form u oc k z . Eq.(17) differs 
from that in classical dynamical critical phenomena in only two ways: the substitution of 
(p — Pc) at T = for (T — T c ) and that the fluctuation frequencies u also has a scale which is 
simply the temperature of measurement T, i.e. the distance along the frequency axis from 
the quantum-transition at (p = p c ) and T = 0. These two ensure that there is a scale T x (jp) 
which gives the crossover from quantum fluctuations at low T to classical fluctuations at high 
T and which marks the temperature below which fermi-liquid properties are expected. As 
opposed to classical critical phenomena, the critically-scaled behavior in physical properties 
is to be expected in the entire regime bounded by the transition temperature T c (p) and the 
crossover scale T x {p) and the upper cut-off temperature of the fluctuations, given by the 
microscopic energy scales (for example, the exchange splitting in a ferro- or antiferro- qcp). 
This is usually referred to as the quantum-critical regime of the phase diagram. 

It is useful to take £ r , & as dimensionless, the former normalized to a lattice constant and 



the latter normalized to the upper cut-off of the fluctuations co c . Eq. (17) has been written 



so that Xo is the spectral weight of the critical fluctuations near the critical point, by which 
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I mean that x"(<2S w,p, T) integrated over all uj and k is xo- The system of-course has other 
fluctuations but xo saturates only part of the total spectral weight. 



The form ( 17) is suitable for discussing experiments as a function of (p— p c ). For analyzing 
experiments as a function of frequency and temperature, an alternate form may be more 
useful: 

*r(q,u>,p,T) =XoT~ d ^(^j- z ,^,^). (20) 

Thermodynamic and transport properties near a Gaussian qcp depend on the cut-off scale 
u c , the spectral function xo> t ne temperature T and the exponent z. The conventional 
theory of quantum criticality of ferromagnets gives z = 3 for ferromagnets and z = 2 for 
antiferromagnets. In this situation, the frequency scale as well as the momentum scale of 
the fluctuations — > as the critical point is approached just as near classical critical points; 
at a temperature T at p = p c , the frequency-scale of the fluctuation is T, and given in terms 



of £t by (20) away from it. 

Let us look at what such critical fluctuations do for T c through what we have just learnt 
in the previous section. In the quantum-critical regime, the scale of fluctuations just above 
T c or the cut-off scale in the fluctuations exchanged is T c itself. This is bad for two reasons: 
the prefactor of the expression for T c goes down with the cut-off, and as discussed and 
illustrated for £ ^ pairing in Fig.Q, inelastic scattering has a particularly bad effect on 
T 

Is this compensated for by increase in the dimensionless coupling constants A^? In other 
words does the peaking of fluctuation at low frequencies and in a region of width around 
the critical wave-vector increase the coupling constants. The answer is sadly, no. The point 
is that the coupling constant depends on the integrated value of the fluctuations over q 
and u and these are essentially fixed by sum-rules. There is a weak (logarithmic) increase 
in the coupling constant if the spectra peaks at low frequencies which is generally quite 
uninteresting. 

Let us consider the role of xo> the spectral weight. It may be taken from the ordered 
moment < M > far away from p c , where it is a slow function of p, to be ~< M > 2 . This is a 
fraction xo of the totals spectral weight of such fluctuations, the rest remaining as incoherent 
particle-hole fluctuations of fermions. As discuss below, the irreducible vertex is reduced by 
Xo compared with the case for phonons. 

In Sec. IV, we will review the firm experimental evidence for the case of the Cuprates and 
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for the Heavy-Fermions that the qcp is not of the Gaussian kind. The fluctuations appear 
to be local in space and decay as power laws in time. In other words x" {q, w,p, T) at p — p c 
has no singularity as a function of q but has a singularity as u and T — > 0. The upper cut-off 
in the fluctuations is at and energy u c , which is only about an order of magnitude smaller 
than T c . Critical fluctuations proposed on phenomenological grounds for the Cuprates have 
the form x" (g, o;,p c , T) = x tanh(w/T), with a cut-off at u c . One could formally put a 



dynamical exponent z — > oo and put this in the form in Eq.(20). But this obscures that quite 
different physical ideas are required to derive such fluctuations than those of the Gaussian 
kind. For example, the derivation of this class of fluctuations ([9], JTU]) for Cuprates rests on 
finding two classes of orthogonal topological variables one of which has spatially local and 
temporally logarithmic interactions while the other has spatially logarithmic interactions 
and temporally local interactions. Classical statistical models are known where the singular 
fluctuations are determined by topological excitations, for example the Kosterlitz-Thouless 
transition in the 2d x-y model and transitions in several vertex models [32] , which do not fall 
under the purview of the classical model of phase transitions which are the inspiration for the 
theory of Gaussian qcp's. The class of microscopic models where such quantum-criticality 
occurs is not known. I will refer to such qcp's as topological qcp's. 

The absence of a spatial scale and the freeing of the scale of frequency of the critical 
fluctuations from the requirement that they tend to lower values as T — > 0, removes the two 
deleterious effect of the Gaussian critical fluctuations on T c : the prefactor in the expression 
for T c remains at u c and there is essentially no extra pair-breaking due to inelastic scattering 
for bj c » T c . The locality of the spatial scale, which in practice means that the spatial 
scale is similar to kp 1 has the additional consequence that the t = and £ ^ couplings are 
similar in magnitude. In fact, given a total spectral weight xo and the I in which pairing 
is favored, it is not possible within Eliashberg theory to have a better spectrum for high T c 
than such a spectrum. 

D. Calculations of T c for the Hubbard Model 

Following the suggestion by Anderson [33] that the essential physics of the Cuprates is 
described by the Hubbard model, there have been innumerable attempts by a variety of 
methods to calculate the properties of the Hubbard model, including T c . The approximate 
methods (RPA [34J, varieties of dynamical mean-field theories [331 l3o]. variational Monte- 
carlo [371 EH]) m comparison with the best numerically precise method, Monte-Carlo without 
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sign problems for U/t up to 7, [39], illustrate the hazards of the enterprise. While the last 
give the upper limit, T c /t < 10~ 3 , the various approximations give a value one and sometimes 
two orders of magnitude larger. 

Since (apart from the RPA), such calculations present results for T c without providing 
the form of the fluctuations spectra or of the coupling function or calculate normal state 
properties, it is hard to say what exactly goes wrong in even such elaborate and careful 
calculations. A hint for one of the possible deficiency from the best Monte-Carlo [39] calcu- 
lations done at various sizes of the lattice (up to 10 x 10) is that the nearest neighbor static 
spin-correlations of the d-wave superconducting wave-functions on a square lattice match 
those of the AFM wavefunction. So, it is possible for calculations on clusters of small size 
to lead to misleading conclusions which disappear for large enough size. 

In a later section, we will review the results of s-wave pairing in the attractive Hubbard 
model, where reliable results without numerical difficulties are available in special limits. 

E. Excitonic and similar Pairing 

Noting the prefactor in the BCS expression for T c , it is tempting to suggest attractive 
interaction among electrons near the fermi-surface through exchange of fluctuations whose 
frequency is high unlike the case of phonons where it is limited by M -1 / 2 , where M is the 
ionic mass [4*0Tf4*2] . (Note that there is no factor of M in the coupling constant A.) These 
other excitations must then be fluctuations of electronic degrees of freedom themselves or 
photons. For the latter, the coupling constant A is inevitably proportional to the fine- 
structure constant. So we need consider only the electronic excitations. It seems that starting 
with Little [10] no known excitation has not been thought of in this context. Although some 
of the ideas are not correct to begin with, most are in principle correct; the problem is in 
the smallness of the coupling constant. When they are wrong in principle, the mistake is of 
the following kind. Consider electron-electron interaction to second order in the Coulomb 
interaction: it is attractive. But if the whole series is summed, it gives just the screened 
repulsive Thomas-Fermi interaction. An example of this in connection with a proposal [H] 
was provided by Inkson and Anderson |43j and further checked through detailed calculations 
[44j . The original proposal by Little [ID] to use the excitons of an insulating side-chain in 
an organic metal has similar problems. 

It should be noted however that essentially every e-e mechanism we can think of is an 
excitation in the particle-hole channel and that if we get a decent T c , it is generally always 
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connected with a high frequency cut-off of the fluctuations. So all e-e processes we consider 
may be called Excitonically induced pairing. But all empirical information on high T c 
induced through e-e processes indicates that the 'excitons' must be particle-hole fluctuations 
of the same one-particle excitations that pair up and such that they engender singularities 
in the single-particle spectra so that the distinction between particle-hole fluctuations of the 
fermions and the 'excitons' is lost. 



III. ELECTRON-PHONON INTERACTION PROMOTED T c 
A. McMillan's Expression for the Coupling constant 

For a general spectral function, 

5 a)/3i7)5 (q,i/,ft) = / du — ^ , (21) 

J-oo il — u 

the Eliashberg equations require a numerical solution. For the limited purposes of this 
paper, it is enough to use the McMillan simplified solution [25], where by the transition 
temperature in the s-wave channel is given for the e-ph problem by 

T c ^cexp(-(l + A e - p,l )/A e - p/l ), (22) 

with 

a- . m<p> t (23) 

where N(0) is the density of states of un-renormalized electrons at the fermi-surface, < u> > 
is an average over the phonon frequencies, k is a (suitably defined) average over the lattice 
stiffness, k = M < u 2 >, and < I 2 > is the scattering averaged over the spectral function 
of the phonons: 

< p >= I dS * I dS * Kql^k ^k' 1 ^ 

J dS k J dSwv^v^ 

where a spin-trace (in the singlet channel) of gr^(k, k + q)g" s (—k, — k — q) has been taken 
to obtain |<?£ q | 2 . 



It is important to note that the A of Eq.(23) also gives the mass renormalization in the 



normal state due to phonons so that the renormalization of the specific heat coefficient due 
to interaction with phonons is 

(m*/m) ph = (l + X e - ph ). (25) 
27 



The coefficient of resistivity due to electron-phonon interactions in the normal state is also 
given in terms of A. For example for temperatures comparable to an larger than u, the 
resistivity due to e-ph interactions in a metal may be written in terms of the scattering rate 

t-\T) = 2n\ e - ph T. (26) 



B. Empirical Relations in transition metal superconductivity 

Study of empirical relations among the parameters determining Transitions temperatures 
in superconductors gives useful insights to the physics of metals. In his analysis of the 
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FIG. 5: Empirically determined "constancy" of iV(0) < I 2 > in classes of metals and compounds. 

properties of superconducting metals and compounds, McMillan [25] noted that iV(0) < I 2 > 
varies only by about factor of 2 while iV(0) and < I 2 > vary by a factor of about 10. However 
the empirically "constant" iV(0) < I 2 > has different values for different classes of metals 
and compounds. For example, see Fig. ([5]), it is is close to one value for the pseudo- 
potential metals like Sn, Pb, Bi and their alloys and another for the transition metals and 
their alloys and yet another in the A-15 compounds. McMillan proffered no explanation 
for the transition metals but showed using the fact that the ion-plasma frequency Q p are 
always much larger than the phonon frequencies, that for the pseudopotential metals that 
iV(0) < I 2 > /Q 2 is approximately constant. 
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Barisic, Labbe and Friedel [15] presented a simple and strong argument for transition 
metals and compounds on the basis of the tight binding representation of the band-structure 
and of the electron-phonon coupling that N(0) < I 2 > is related simply to the cohesive 
energy E c of the metal. The argument is summarized in appendix A with the conclusion 
that 

N(0) < I 2 >^ N(0) < d 2 E c /dR 2 >^ N(0)E c /r 2 . (27) 

< d 2 E c / dR 2 > is the average of the second derivative of the change in kinetic energy of the 
metal as the nearest neighbor distance between two atoms R is changed leaving the others 
fixed; ?"o is the size of the typical metallic orbital. 
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FIG. 6: Empirical Relation between experimentally deduced A and the bare electronic density of 
states at the chemical potential in 3d-4d and 5d metals and their alloys. 

One may be tempted to conclude that since N(0) < I 2 > within a given class of transition 
metals or compounds is approximately a constant, one may simply increase A by reducing the 
average lattice stiffness M < u> 2 > and thus increase T c . This led to the soft-phonon myth, 
much propagated in the 1970's . Quite apart from the fact that the prefactor < u > would 
prefer matters the other way for high T c , there is also another empirical rule [16J H?j followed. 
It is that within a given class of materials < I 2 > f(M < u 2 >) is also approximately a 
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constant. For transition metals and alloys this is exhibited in fig. ([6]). Only the first half 
of the 3d metals is superconducting; for 4d and 5d metals the two lines correspond to the 
fermi-level in the bonding part of the band and in the anti-bonding part of the band. This 
"constancy" with a different value is also followed in the A-15 family of superconductors 
[17] . This rule may also be derived within the tight-binding approximation [15] and is 
summarized in the appendix. The final result is that for N(0)W » 1, where W is the 
electronic bandwidth 



This is in qualitative accord with the data in Fig. ([6]); the 3d, 4d and 5d metals have 
progressively larger band-widths and the anti-bonding part of the band has effectively a 
larger band- width than the bonding part. 

1. Maximum T c and General Lessons. 

It is an amusing excercise to use the empirical relations to see the scale of maximum 
T c . Such an excercise was indulged by McMillan |25j and by Anderson and Cohen [16] for 
pseudopotential metals with amusing misunderstandings resulting, for which the originators 
can hardly be blamed. The emphasis in these estimates was both on conditions on lattice 
stability and the compeition between the repulsive Coulomb pseudopotential and the at- 
tractive interaction through phonons. I give here an estimate for maximum T c based on the 
empirical relations discussed above, which is similar to that done about 30 years ago [T5] . 
This is valid for metals and compounds where a tight-binding approach to the electronic 
structure and electron-vibrational interaction is valid. These are the likely materials for the 
highest T c by e-ph. 

Using the expression, T c ~< to > exp —(1 + A)/A, and assuming < cu >^< cu 2 > 1 / 2 , one 
can use the two empirical relations to express T c in terms of A, E c /r^ and the ion mass M. 
Now maximizing T c with respect to A gives that the optimal value of A ~ 2 and 



One may estimate the magnitude for compounds with Nb as the main element with atomic 
mass 92. The typical cohesive energy per formula unit estimated from the second moment 



A = N(0) < I 2 > I (M < to 2 >) oc 1 + N(0)W 



(28) 




(29) 
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of the band is about 10 eV for Nb or Nb^Sn. Taking tq = 1.34A, the co-valent radius of Nb 
gives (T c ) max « 30K. This is to be compared with the experimental value of 9.2 K for Nb 
and 22 K for Nb 3 Sn. 

The conclusion is that for the highest T c from electron-phonon interactions, one needs 
a small mass M of the ion to give the large frequency scale of attractive interactions. In 
a compound, the electrons near Ef must have substantial weight on the same atoms with 
the small M so that the electron-phonon interactions and the phonon energies are related 
to the average stiffness Wjr\. The highest T c of the electron-phonon variety so far is MgB 2 
which satisfies these conditions. Taking the geometric mean of the mass ~ 13 proton mass 
and W from band-structure of about lOeV and r ~ lA gives a (T c ) max of about 60 Kelvin 
while the actual value is 42 kelvin. One cannot do better than metallic Hydrogen where the 
bandwidth is estimated to be about leV. (T c ) max of about 100A" is to be expected. These 
are generous limits because they exclude the reduction due to the Coulomb pseudo-potential 
which is expected to be especially severe for hydrogen. 

An important lesson from the study of the empirical relations in "high-temperature" 
superconductors of the e-ph variety and their explanations, is that the parameters deter- 
mining T c are gross parameters, which depend on the average local stiffness of the lattice 
and on how the electronic energy of the bands changes with local deformations. Although 
superconductivity is a fermi-surface phenomena, the parameters determining T c are prop- 
erties related to the variation of the electronic bonding energy with the local fluctuation 
responsible for pairing. These lessons carry over to the electronic mechanisms of pairing 
except that getting effective coupling constants of 0(1) with their main benefit - the larger 
high-frequency cut-offs, requires rather more stringent conditions, as we shall see below. 

IV. SUPERCONDUCTIVITY FROM FERMION INTERACTIONS 

The theory of Fermi-liquids by Landau [48] . which was almost concurrent with the BCS 
theory of superconductivity [TJ, led to an enormous interest in the experimental study of the 
properties of liquid He 3 in the 1960's [45] . Liquid He 3 near the melting line is very strongly 
correlated with the magnetic susceptibility about 25 times larger than a non-interacting 
fermi-gas of the same fermi-energy. With BCS theory in mind, it was natural to think of 
pairing or superfluidity in He 3 . Two He 3 atoms have a large s-wave repulsion due to the 
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hard-core interaction. It is necessary therefore that any bound state of a pair of atoms ip{ r ) 
have a node in the wave-function at the relative co-rodinate r = 0. The radial part of the 
pair wave-function must therefore vanish as ip(r) oc r e , and so the bound-state must be in 
the £-th angular momentum. It follows from Pauli-principle that even £'s have total spin 
zero, i.e. a singlet state, while the odd ts have spin 1, i.e. a triplet state. The idea of a 
finite angular momentum pairing for He 3 therefore arose [50] . 

Following Berk and Scrieffer's [51] result that ferromagnetic fluctuations are deleterious 
for s-wave pairing, Layzer and Fay [52] suggested that such fluctuations may promote spin- 
triplet I = 1 pairing. The discovery of such pairing in liquid He 3 , (for reviews, see 13], 153]) . led 
to further theoretical ideas and calculations. It is generally agreed that although the idea of 
exchange fluctuations give the correct symmetry of pairing, weak-coupling calculations based 
on the idea do not work quantitatively. Calculations properly taking into account the short 
range repulsion between the He atoms and estimating interactions with constraints put by 
the measured Landau parameters give the right scale of T c and its pressure dependence [54J . 

In the late 1970's superconductivity in heavy fermions was discovered [551157] and more 
superconducting heavy fermions continue being discovered. In these materials, the mass 
enhancement of the fermions is of O(10 3 ), so that the effective fermi-energy is of the same 
order or smaller than the characteristic phonon energy. It was suggested that in this case 
the phonon attractions could not overcome the Coulomb repulsion because the concept of 
the Coulomb pseudo-potential is invalid and therefore the pairing must be in a finite angular 
momentum state [58] ■ Experiments were suggested to test this suggestion which were soon 
carried out [59|. Transport experiments were analyzed j6(Jl El] and they could only be 
understood if there was a line of zeroes of the gap-function on the fermi-surface. Such a 
gap function is not allowed [62j [63] in the spin-triplet manifold in the presence of spin-orbit 
scattering. Therefore i = 2 pairing was to be expected. At the same time, it was known 
that the heavy fermion superconductors were close to anti-ferromagnetic instabilities. This 
led to an investigation of pairing through anti-ferromagnetic fluctuations which we have 
reviewed above [T9] which has since been much used (and in my view abused) extensively 
in connection with the superconductivity in the Cuprates. A concurrent RPA calculation 
[61] of the pairing susceptiblity in the Hubbard model also revealed a tendency to d-wave 
pairing. 
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A. Superconductivity in the dilute Fermion Gas with varying Interactions: 
Theory and Experiments for Cold Atoms 



The advent of the technique of cooling atoms in atomic traps has generated (among other 
things) a new class of Fermi-liquid in which the particle density v is low but the inter-particle 
interactions U can be varied from very weak to very strong. The effective particle-particle 
interaction in the low density limit is completely specified by the scattering length a so that 
all physical properties are functions of the dimensionless coupling strength 

k=^-, (30) 
k f a 

where kf is the magnitude of the fermi- vector. The theoretical deduction of the effective 
interaction in this situation is simpler than for other fermi-liquids. The upper cut-off of 
interaction energies is the Fermi-energy. The essentially exact calculations of T c /Ef possible 
for this problem provide a measure of the highest value to be expected from e-induced s- 
wave pairing for the complicated situations. It is not coincidental that the maximum T c /6d 
realized in e-ph superconductors approaches the highest value in the calculations summarized 
below. 

The t-matrix approximation for effective interactions, which is exact in the low density 
limit, gives that the the scattering length a is given in terms of the parameters of the 
Hubbard model by [5] 



Here U is the interaction parameter in the Hubbard model, < to give s-wave pairing, 
e(k) = k 2 /2m, and U~ = —Tilom; 2ti/Iq is the upper cut-off in the integral over k introduced 
to avoid the ultraviolet divergence. 

The weak-coupling limit (BCS-limit), when the attractive interaction strength is negligi- 
ble compared to the kinetic energy is given by k — > — oo and the opposite or Bose molecular 
limit by k — y +oo. In between is the unitarity limit k — > where a — > oo. Universal results 
are to be expected for physical properties for all low density attractive interaction models 
in these three limits. These limits are realized in experiments by tuning the interactions 
through the Feshbach resonance. 

Monte-carlo methods have been used to calculate [5] T c /Ef both for interactions for free 
fermions around the unitarity limit as well as for attractive interactions in the Hubbard 




(31) 
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FIG. 7: Monte-Carlo calculations [5j for dilute fermion gas with varying interaction parametrized 
by the dimensionless number k around the unitarity limit, n = . The results are compared to 
the weak interactions or BCS limit, k — > — oo and to the Bose molecular formation limit, k — > oo. 
Non-universal results are obtained for the latter limit, shown by different values for different forms 
of potentials 

model on the square lattice for varying density. Figure ^ gives the results and compares 
them to the results in the BCS limit as well as the Bose molecular limit. The results in 
the unitarity limit are T c /Ej = 0.152(7). The same result is obtained in the unitarity 
limit for the attractive Hubbard model, as expected from universality. Figure ([7]) shows, 
unexpectedly, that the results as a function of k are not monotonic in going from the BCS 
to the molecular limit, but have a maxima on the molecular side. 

The thermodynamic deductions of the properties of the cold atomic gases in the unitarity 
limit [6] give a value T c /Ef ~ 0.2. Given the difficulty in determining thermodynamic 
properties precisely as well as the non-uniformity in the density of the gas in the optical 
trap, this should be considered good agreement with the calculations. 

B. Liquid He 3 : Pairing Symmetry, T c and Connection to Landau Theory 

Near the melting line of He 3 , the effective mass is about 6 times larger, the compressibility 
about 15 times smaller and the magnetic susceptibility about 25 times larger than a non- 
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interacting gas with the same density of states at the chemical potential. The transition 
temperature of superffuid He 3 might therefore be taken to represent the scale of T c to 
be expected for a strongly interacting fermi-liquid away from criticality. The low energy 
properties are given by the Landau theory and thermodynamic and transport properties 
have been measured extensively. T c /Ej is about 10~ 3 . Much is to be learnt from the work 
on this problem, not the least is why T c is so low. Variation of T c with pressure together 
with results of an interesting calculation by Patton and Zarnalingham (65] are shown in fig. 

The general pairing vertex has been given in Eq. (|2| and Fig. ^ in terms of the irreducible 
particle-particle vertex 1$. The relation of Is(k, k + q, k') to the Landau parameters is in- 
teresting to explore. With k, — k and k + q, k — q, all close to the fermi-surface, is(k, k + q) 
is specified by Is(&, 4>), where 9 is the angle between kp and k + q F and <fr is the the angle 
between the planes formed by k, k + q and k', k' — q. For pairing, we need only Ig(7r, 0). 
The forward scattering limit q — > 0, u — > of I are the province of the Landau theory of 
Fermi-liquids. The "A" Landau interaction parameters are given by the forward scattering 
limit, (vfq)/u — > 0, then oj — > of N(0)Is(9, 0) = As(tt,4>)- They are decomposed into 
different angular momentum channels: 

4 = A 5 (vr,0)P,(0). (32) 

The forward scattering sum rule (imposed by Pauli-principle) requires 

J2 A s,e = °- ( 33 ) 

s,i 

Patton and Zarnalingham [65J argued that since the irreducible pairing vertex needs to be 
calculated in the limit cu/lvpq) « 1, it should be related to the "F" Landau parameters, 
which are in related to the irreducible vertex in the other limit, u/(vFq) —> and then 
q -)■ 0. "F's" are related to the "A's" through 

F s>i = A S;i /{l + A s/ ) (34) 

The measured compressibility and the spin-susceptibility at various pressures provide F 0i o 
and Fi t Q respectively as a function of pressure, while the specific heat provides Fq^. The 
Landau parameters are expected to rapidly decay with increasing i due to the short-range 
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FIG. 8: T c in liquid He 3 as a function of Pressure compared with the calculations in Ref. (|65|) based 
on fitting Landau Parameters 

nature of the effective interactions; if one assumes that it is saturated by £ = and £ = 1, 
Fi t i may be extracted using the forward scattering sum-rule. 

Patton and Zarnalingham [65J found that the pairing interactions extracted with the 
assumptions above give repulsion in the £ = 0, s = and attraction in the £ = l,s = 1 
channel. They could get a surprisingly good value for T C (P) over the whole range of pressures 
(see fig. ([8])) using the BCS formula T c m w c exp (— 1/A) if they use that the upper cut-off 
"c(P) « ±E f {P). 

The weak part in the argument is the assumption that uj/Vjq « 1 may be used with 
q — > to related the vertex to the "F's" and the lack of any explanation of the value of 
the cut-off. Detailed calculations of the Landau parameters and the extension of Landau 
theory to vertices with finite momentum transfer (to remove one of the the weak parts of 
the argument) have since been done using partly phenomenological and partly microscopic 
approaches [3]. The complete such calculations done by Pfitzner and Woelfle [51] carefully 
using the exact equations for the vertices and using empirical information on the measured 
thermodynamic and transport properties of the normal state of liquid He 3 to constrain 
limiting values of the solutions and interpolating to get the complete vertex. This is the best 
kind of microscopic theory suitable for such difficult problems. With these calculations the 
singlet £ = 0,2 channels are shown to be repulsive and the triplet £ = 1 channel is attractive 
with the right variation of coupling constants with pressure to provide the variation of T c . 
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One of the important lessons to be learnt from such studies is that the effective coupling 
constants are determined by the A(_ parameters (and their finite momentum extensions which 
are expected to be smoothly varying) and their modulus is always smaller than 1. This 
represents the physical fact that the dimensionless coupling constants involve the product 
of the density of states and the interactions. Large renormalizations always (nearly) cancel 
out in this product because of cancellations of self-energy and vertex renormalizations. In 
Landau theory these cancellations are related to conservation laws. The other lesson is the 
empirical lesson that u c is an order of magnitude smaller than Ef. This cannot be obtained 
from Landau theory but any reasonable microscopic calculation gives that the cut-off of the 
fluctuations goes down as the F parameters increase. For example the sound-velocity goes 
down as the compressibility or Fq increases, the spin-fluctuation frequencies go down as the 
magnetic susceptibility or (m* /m) (1/1 + Fq) increases, etc. 

Let us turn briefly to the microscopic calculations from weak-coupling theory. In the 
paramagnon model, the best of these is due to Levin and Vails [21]. The local Hubbard 
interaction 



promotes ferromagnetic exchange in free-fermions (no lattice potential) and associated in- 
crease of the amplitude of ferromagnetic spin-fluctuations for / < I c , the critical value for 
ferromagnetism. This model does not describe the physics of He 3 very well because of the 
substantial range of the hard-core interactions compared to inter-particle spacing; it does not 
give the right values of the Landau parameters or their pressure dependence. But neverthe- 
less several features of our interest in relation to calculations of T c , from fluctuations induced 
by particle-hole interactions, from such a model (and its modifications) have substantial ed- 
ucational value because RPA respects conservation laws. Levin and Vails performed such 
calculations in the Hubbard model to calculate both the Landau parameters as well as solve 
the Eliashberg equations with the fluctuation spectra obtained in the calculations. The 
important results of their analysis and calculations are: 

(1) The pairing interaction in the t = 1 channel V\ has a direct dependence on I but so 
does the effective mass as well as a m*/m. The renormalized pairing interaction parameter 
Aidepends on the product V\(m* / m) _1 . Moreover the self-energy correction in the T c formula 
depends on parameter Ao- This goes up with /. This reduces T c and may be taken to 
contribute to an effective reduction of cut-off if one insists on using the BCS formula for T c 




(35) 
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even for finite I. 

(2) The net effect still is that T c goes up with increasing / except close to the ferromagnetic 
transition, where it swings downwards towards 0. This is due to the pile up of the fluctuation 
spectra to low energies for / close to I c . This has two deleterious effects on T c : it increases 
inelastic scattering and reduces the cut-off u c . 

These conclusions are consistent with those from the analysis of the Eliashberg equations 
with a general form of pairing and self-energy kernel discussed above as well as with the 
Ward identities. A Gutzwiller-type variational wave-function with a Hubbard model on a 
lattice j3] with less than 1/2 filling gives much better account of the ground state properties 
of normal He 3 than does a paramagnon model. But the extensions to the excitation spectra 
with the same basic physics can not work, since putting the model on a lattice promotes 
antiferomagnetic correlations and does not lead to pairing interactions in the £ — 1, triplet 
channel. 

C. Experimental results on Heavy Fermions 

Heavy Fermion compounds show superconductivity, generally associated with an AFM 
quantum-critical point, (but note the interesting case of CeCu^Si^ under pressure, where 
two forms of criticality, AFM, and Mixed- Valence each seem to have an associated supercon- 
ducting region [66]), although the converse is not true; AFM quantum critical points in some 
heavy-fermion compounds are not accomapanied by superconductivity. Superconductivity 
always appears to be in a finite angular momentum state and is not due to electron-phonon 
interactions. 

Heavy Fermions in rare-earth and actinide compounds are the ultimate realization of the 
ideas of analyticity and continuity which underlie the Landau quasi-particle idea. In their 
fermi-liquid regime, the effective mass enhancement in several heavy-fermion compounds is 
0(10 3 ) and the quasi-particle renormalization residue z is 0(1O -3 ). This situation changes 
in the quantum-critical regime where the quasi-particle idea breaks down and transport and 
thermodynamic properties are not those of a Fermi liquid. This is a beautiful example of 
how z — > as T — > 0, only logarithmically, produces completely different physical properties 
at low temperatures than z ~ 10 -3 . 

Knowing the fermi-liquid renormalizations is not as useful to deduce parameters for su- 
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FIG. 9: Phase Diagram of a 115 compound and alloys at various pressures shown. Figure is 
reproduced from Ref . (j67|) . 

perconductivity in heavy-fermions as in liquid He 3 for two reasons: superconductivity is near 
qcp's, where such renormalizations are singular and the qcp are at large g-vectors, where 
Landau parameters are not defined. However, the energy scales are set by the renormaliza- 
tions given by the Landau parameters above and are therefore essential to bear in mind. 
They are ideal systems to study magnetic fluctuations by inelastic neutron scattering. But 
only in a few of them are such results available near quantum-criticality because for the 
technique to be fully effective requires large single-crystals. From the study of the ther- 
modynamic and transport properties (such as residual resistivity, temperature dependence 
of resistivity, ultrasonic attenuation, thermal conductivity and nuclear relaxation rates) in 
the fermi-liquid regime, the leading renormalizations in heavy fermions were found to be 
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FIG. 10: T c vs. lb, a measure of the Fermi-energy obtained from magnetic susceptibility measure- 
ments from several heavy-fermion superconductors as well as Cuprates; figure is reproduced from 
Ref.(l67|). 

qualitatively different from that in liquid He 3 . The renormalizations are characteristic of a 
sub-set of fermi-liquids in which the single-paricle self-energy is very weakly dependent on 
momentum compared to on energy. This idea gives that [68J 

m*/ m = 1/z = (1 + F S ) (36) 
K / K ° = ( m */m)/(l + F S ) = 0(1) 
X /xo = K/m)/(l + F a ) = 0(l). 

The enhancement of the specific heat of O(10 3 ) gives z of O(10 -3 ); similar enhancement 
of susceptibility gives Fq of 0(1); the lack of renormalizations in ultrasonic attenuation 
and zero-temperature resistivity give that Fq w 1/z. The temperature dependence of the 
resistivity in the Fermi-liquid regime is oc T 2 with its coefficient renormalized by 0(1/ z 2 ) 
which also follows from the momentum independence of the self-energy. These ideas are also 
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realized in microscopic calculations based on dynamical mean-field theory [691 [70], which 
also start with the assumption that the self-energy is momentum independent. 

At a microscopic level, it is understood that the renormalizations in the heavy-fermion 
lattice are closely connected with the Kondo effect of the magnetic ions in the bath of 
conduction electrons. This physics is well-incorporated in the models of the lattice by 
a variety of different techniques, the most versatile of which is the dynamical mean-field 



theory which yields results similar to Eqs. (36) together with the frequency dependence of 
the one-particle spectra. A theory of two interacting Kondo impurities [7T] with interactions 
among the magnetic ions competing with the Kondo effect at either also exists. The theory 
yields a low energy effective Hamiltonian for the problem with several interesting terms, 
including pairing interactions among the quasi-particles. A low energy Hamiltonian for the 
lattice at a similar level has not yet been derived but is essential for a detailed theory of 
superconductivity in them. As discussed above, one can motivate the d-wave symmetry 
of superconductivity in most of the heavy-fermions through exchange of antiferromagnetic 
fluctuations. But this does not answer the question of the energy scales of fluctuations in 



the critical regime. All we know is just that the renormalizations discussed through Eqs. (36) 
give the scale of the upper limit of the cut-off energy of the pairing fluctuations to be down 
by z from the bare parameters, i.e. of O(Ef), the renormalized fermi-energy. (As discussed 
already, the value of z serves to renormalize the cutoff scale as well as the coupling constants 
in general accord with Landau fermi-liquid theory.) The actual experimental values of T c 
are down from this by only between one and two orders of magnitude from the upper cut-off, 



see Fig. (10). One may speculate that this is associated with the fact that unlike liquid He , 
superconductivity in heavy-fermions generally always occurs in the region around quantum- 
criticality, see Fig. ([9]), and the criticality is not of the conventional variety where the energy 
scale of fluctuations — > as criticality is approached. 

In the one case in which quantum-critical fluctuations have been studied in single-crystals 
by neutron scattering the criticality is of the topological or local variety [72] discussed in 
Sec. II: the fluctuation spectra scales as u/T but is essentially q- independent; the space and 
time are not tied together as in Gaussian quantum-criticality. This particular compound 
is however not superconducting. But, based on the measured specific heat, and resistivity 
measurements, one can be fairly sure that criticality is of the same variety in the heavy 
fermions where superconductivity occurs with relatively large values of T c j Ef. The measured 
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specific heat oc TlogT and temperature dependence of resistivity oc T cannot be understood 
in three-dimensional materials by the Gaussian variety of quantum criticality [3T] ■ 

The attempts to derive quantum-criticality in heavy fermion lattice [73] are equivalent to 
the self-consistent solution of the single-impurity Kondo problem in a conduction electron 
bath which is close to antiferromagnetic quantum-criticality [71]. The breakdown of the 
Kondo effect near magnetic singularities is doubtless an important aspect of the problem but 
such methods do not represent a solution of the criticality problem. The issue of quantum- 
criticality in heavy fermions is a beautiful largely unsolved problem; it is also not possible 
to discuss their "high" T c without understanding the fluctuation spectra near quantum 
criticality and its coupling to fermions. At a phenomenological level also, one is limited in 
this discussion due to lack of measurements of the fluctuation spectrum or single-particle 
spectra in the heavy-fermions near quantum-criticality. Such measurements will require 
quality and size of samples similar to some of the Cuprates. 

D. Superconductivity in the Cuprates 




FIG. 11: The Universal Phase diagram for hole-doped cuprates based on properties which show 
characteristic changes across the lines drawn in all Cuprates. The demarcation region between 
Region II and the AFM region has not been determined clearly by experiments or theory. 

The family of Cuprates have the highest transition temperatures discovered so far. 
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This phase diagram of Fig. (11) organizes the extensive experimental investigations in the 
Cuprates in the last twenty years. It is the universal phase diagram in the sense that the 
lines drawn demarcate change in properties in the same manner in every Cuprate family 
investigated. (There are of-course Cuprates with special properties in one or the other re- 
gions of the phase diagram which are not shared by other Cuprates.) The phase-diagram is 
meant to be schematic and correct only in its topology. 

There have been innumerable ideas and calculations to understand the properties of the 
Cuprates. I will only summarize the work which seeks to understand all the principal normal 
state properties as well as superconductivity with one set of ideas and has predicted unique 
results which have been tested by experiments. 

The central organizing feature of the phase diagram is a quantum-critical point from 
which two lines emanate which separate the normal state in to three parts; The region of 
superconductivity is located around the qcp. The location of the qcp within the diagram 
varies with respect to the region of T c ^ in different cuprates; the maximum T c is generally 
not at the doping x c of the qcp. In Region I, quantum-critical properties are observed in Re- 
sistivity, optical conductivity, Raman Scattering, nuclear relaxation rate on Cu-nuclei (but 
not O- nuclei which show fermi-liquid type relaxation rates), tunneling into a conventional 
metal and single-particle spectra measured by angle-resolved photoemission (ARPES). The 
superconductive pairing is in the "d-wave" channel. The great mystery has been the nature 
of Region I and of Region II. Region II, the so-call pseudogap region abuts the bound- 
ary to Region I whose properties are determined by quantum-critical fluctuations (qcf). 
Through extensive experimental investigations, it is now generally believed that Region II 
represents a phase which competes with superconductivity, much the same as antiferromag- 
netism competes with superconductivity in the heavy fermions. Needless to say, Region II 
has no long-range anti-ferromagnetic order and near the region of maximum T c , in at least 
one of the heavily investigated Cuprates, the magnetic correlation length is about a lattice 
constant, which is the scale of magnetic correlations in non-interacting fermions. 

Given the general proposition that the pairing symmetry and T c are the properties of the 
normal state just above T c , the central problem for superconductivity in the Cuprates is the 
understanding of the nature of the quantum-critical fluctuations, which is a problem which 
can be usefully attempted only if the nature of Region II is understood since given the phase 
diagram, the fluctuations in Region I are the quantum-fluctuations of the order parameter 
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in Region II. There is physics in both regions which we had not come across before. For 
example, as we will see, the fluctuations in Region I are quite unlike the Gaussian fluctuations 
discussed in Sec. II C. The symmetry breaking in Region II remains controversial. 

There have been two valid reasons for skepticism that Region II has a distinct order 
parameter. One is that it the extensive investigations on the Cuprates with every available 
tool and some newly invented had not discovered any credible broken symmetry and the 
second that although the specific heat decreases near below the line demarcating Region I 
and II, there is no evidence for any singularity in it as in the typical classical phase transition. 
However, in the last five years, an elusive proposed order [75H77] has indeed been discovered 
experimentally in four distinct families of Cuprates [7HH8T] . Moreover the model representing 
the observed classical phase transition to this order has a non- diverging specific heat at the 
transition [32], whose quantitative value is consistent with the experiments [82J. 

I will now summarize the theoretical results briefly. In a two-dimensional Cuprate model 
with three degrees of freedom per unit-cell, representing the dj.2_ y 2-orbital of Copper and 
the p x and the p y -orbitals on the Oxygen atoms on the x-direction and the y-direction 
respectively of the Copper, together with on-site and nearest neighbor interactions, two 
phases which are odd under time-reversal due to spontaneously generated orbital loop- 
currents but preserving the translational symmetry were shown to be locally stable in a 
mean-field calculation. The phase observed has the symmetry of one of these phases but 
there are some other details [78TI8T] which are different and obtainable only in a more 
complicated model [83J . The observed as well as the predicted phase has a symmetry which 
can be characterized as the ordering of a polar time-reversal odd vector L. L describes the 



loop-current pattern shown in Fig. (12) which preserves even-ness only in reflection on one 
of the four-possible symmetry planes of the square lattice. L has four-possible orientations, 
corresponding to the four possible domains of the ordered phase. In the quantum-fluctuation 
regime, the instantaneous pattern of currents is described by vectors Lj, which vary in 
direction among the four-possible orientations. The quantum-fluctuations are described 
by the correlation function [9] [10] < XJf(t)XJj(t') >, where Uj is the generator of the 
rotation among the four-configurations of Lj. It turns out that in the fluctuation regime 
the discreteness of Lj can be replaced by a continuous vector so that the fluctuations are 
related to the fluctuations of a quantum-fluctuations of a 2— dimensional rotor. In a model 
including dissipation due to the coupling of the fluctuations of Uj to the fluctuations of the 
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FIG. 12: The four possible domains with the symmetry of the observed order in the underdoped 
Cuprates. 

fermion-current, the spectral-function x"(q, of the fluctuations in the quantum-critical 
regime has been derived ([9], [ID]) to be of the form 

-Xotanh(w/2T), |w| < u c ; 
bnx(q,w) = ^ (37) 

This is precisely of the form which was suggested [51] to account for the singular transport 
properties in Region I of the phase diagram and is the basis of the marginal fermi-liquid. 
The single-particle scattering rate from such fluctuations can be calculated, the imaginary 
part of the self-energy for frequencies much larger than the temperature is given simply by 




Im£(u;,k) = --A(k) <[ '~ (38) 



Here A(k) is a coupling function whose derivation is discussed below. Fig. (13) shows 
that in the (7r,7r)— direction, this prediction is fulfilled in all Cuprates in which ARPES 
measurements have been carried out. The scattering rate at any energy u is essentially 
oc ^ q _/ x"(q, k/)- Therefore the linearity in uj of the scattering rate up to some energy 



and constancy thereafter is a direct proof of the fluctuation spectra given by Eq. (37). Earlier 



experiments with greater resolution at lower energies providing evidence the crossover from 
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linearity in T to linearity in u have been reviewed [85]. Fig. ( 13 ) is also a proof that a distinct 



fluctuation spectra with a sharp cut-off u c ~ 0.5eV exists universally in the cuprates. These 
critical fluctuations themselves for q — > are directly observed in Raman scattering (but 



the experiments have not been carried out all the way to the cut-off energy), see Fig. (14), 
where evidence for the universality is presented through S(u) = (1 + n(u/T))x" (w) in the 
limit q — > 0. 



The spectra of Eq.(37) is quite unlike the Gaussian critical spectra discussed in Sec.IIC. 



The singularity at (u, T) — > does not affect the bulk of the spectra at all which extends at all 
T to u) c , which as we will infer from experiments is about Ef/4. Second and most curiously, 
the critical spectra has no spatial scale, the concept of a dynamical critical exponent z is 
lost. 
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FIG. 13: The linewidths of the Momentum Distribution curves for different cuprates as a function 
of the energy. The imaginary part of the self-energy is obtained by multiplying this linewidth with 
the bare fermi- velocity. The detailed references for each cuprate are given in Ref. (|86p 

Given that these singular fluctuations determine the properties above T c including the 
scattering rate of the fermions, it is natural to ask if they promote superconductive pairing 
with the observed d-wave symmetry and with the right order of magnitude of T c . As has 
already been noted in Sec. II, the spectra of the fluctuations (eq:flucspec) is ideal for high 
T c on the basis of Eliashberg theory. It has a high upper cut-off, and it has the least 
inelastic scattering possible in a quantum-critical spectra. However, the q-independence 
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of the spectra makes one wonder if it can promote d-wave pairing. To investigate this, 
the momentum dependence of the coupling of the fermions to the fluctuations has been 
calculated. 

The coupling of the fermions to the fluctuations has been calculated [HZj- In the contin- 
uum limit, U(r) is the angular momentum operator generating rotations of L(r). Therefore 
it can only couple to the local angular momentum of fermions. So the coupling is of the 
form 

H mt ex J dr ^# + (r )( 7)(f x p)^(r, a)U(r) + E.G. (39) 

E int has also been derived for the fermions in a two-dimensional model of the Cuprate lattice 
and the coupling constant g estimated in terms of the same microscopic model which gives 
the symmetry of the observed order and its approximate magnitude. It is useful to note 



that Eq. (39) is the natural orbital angular momentum analog of the familiar collective 



spin-fluctuation coupling to to spin-flip excitations of fermions. 



We may write Eq. ( 39 ) in momentum space 



H int = ^ x k')^ + (k, a)V + (k', a)U(k - k') + E.G. (40) 

k,k',o- 

The coupling constant of the scattering of fermions to the fluctuation spectrum can be ex- 



tracted from the ARPES data in the normal state, Fig. (13) [86J. From such measurements, 
one deduces that the coupling constant A for all Cuprates measured by ARPES is between 
about 0.7and 1 and the cut-off u c is between 0.4 eV and 0.5 eV. The lattice generalization 



of Eq. (40) also predicts that the scattering rate varies oc a + 6cos(4#) where 9 is measured 
from the tt, it direction with a variation of about a factor of 2 going from the n, it to the n, 
direction. 

The momentum-dependence of coupling, even though the spectrum itself is momentum 
dependent is crucial to the symmetry of superconductivity promoted by the critical fluctu- 
ations. This is seen as follows: Integrating over the fluctuations in Eq.(40)gives an effective 
vertex for scattering of fermion-pairs: 

E pairing n J2 A(k,k / )i(-k , )4(k , )c CT (k)c CT ,(-k); 



kcrk' 



A(k, k') =1 (k, k'M-k, -k')Re X (u = £k - 4)- ( 41 ) 
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FIG. 14: Universality of Raman scattering in Cuprates near and Ba x K\- x BiOz for compositions 
of doping with the highest T c . The figure is taken from 



This is exact to O(A^y), where A's are the dimensionless coupling constants exhibited below. 
In the continuum approximation for fermions near the fermi-energy, 7(k, k') oc i(k x k'). 
The pairing vertex is then 

A (k, k, ) oc -(k x k') 2 Rex(k - k', u). (42) 

Since Re%(k — < for — uj c < u < u c , independent of momentum, the pairing 

symmetry is given simply by expressing (k x k') 2 in separable form : 

(k x k') 2 = 1/2 [(kl + kl)(k'* + k'*) - (kl - - k'Z) 

- 4(k x k y )(k' x k')] ■ (43) 



Pairing interaction in the s-wave channel is repulsive, that in the two d-wave channels is 
equally attractive, and in the odd-parity channels is zero. The factor i in 7(k, k'), present be- 
cause the coupling is to fluctuations of time-reversal odd operators, is crucial in determining 
the sign of the interactions of the pairing vertex. 

To estimate T c , we use what has been discussed in Sec. (II) about the effect of inelastic 
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scattering on finite angular momentum pairing. T c is given approximately by 



T c «w c exp(-(l+|A a |)/|A < ,|), (44) 

where X s is the coupling constant which appears in the normal self-energy and Ad, the cou- 
pling constant which appears in the pairing self-energy. From the measurements summarized 



above and Eq.(42), Xd/X s ~ 1/2- Using the deduced value of X s and u c from the ARPES 
measurements, one estimates a value of T c »s 100K. 

Although T c is expected to reduce in the underdoped region due to the competing phase 
and in the overdoped region due to the change in the spectra to an incoherent spectra below 
a cross-over scale, no quantitative calculations for these effects exist. 

E. The case of the Fe-Pnictides 

The newly discovered superconductivity in the Fe-pnictides is also quite unlikely to be 
induced through interaction with lattice vibrations. A recent review is Ref. (l89p . The highest 
T c in this class of compounds so far is about 50K in RF eAs(Oi_ x F x ) : (R = Ce, Pr, Sm, Nd, 
etc). The "parent compound" at x = is metallic but antiferromagnetic. The lack of 
significant observable feature in the specific heat at the high superconducting transition 
temperatures raises doubt as to whether bulk superconductivity in this structure of the Fe- 
Pnictides has indeed been found. There is also some evidence that this structure may have 
a two-phase co-existence as a function of doping. A closely related new structure of Fe- 
Pnictides called 122 appears to form good single crystals. Thermodynamic data indicates 
bulk superconductivity. The phase diagram of this class of materials appears similar for 
both hole doping and electron doping. The superconducting region is organized around a 



quantum critical point, see Fig. (15) where an AFM/structural transition temperature — > 
with change in doping. The mystery of the Cuprates: the nature of the ordered phase on one 
side of the critical point is absent. Moreover, the anisotropy in resistivity of these compound 
is less than an order of magnitude - they are properly considered three-dimensional. The 
fermi-surface has five sheets, most prominently a pair of electron-pockets centered at the 
zone-center and a pair of electron pockets at the zone-faces. At this point, the symmetry 
of superconductivity is not unambiguously known. ARPES experiments indicate that there 
is a gap everywhere on the Fermi-surface. This may well be an extended s-wave form of 
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FIG. 15: The Phase diagram of Co-doped FeAs2 [90J, [91J. The green triangles mark the transition 
to an altered structure while the black circles mark the antiferromagnetic transition. Supercon- 
ducting region is shown in blue. 

pairing, as has been suggested [92--94J. As discussed earlier, AFM fluctuations would favor 
such a symmetry of pairing for an appropriate fermi-surface; the form of the fermi-surface 
in this class of compounds is the right kind. Given the phase diagram, it is also reasonable 
to infer that the AFM/structural quantum-criticality is responsible for the high transition 
temperatures. But as discussed above, the Gaussian quantum-criticality is bad for T c , both 
due to the prefactor and the inelastic scattering, especially for the case of pairing not of 
the simple s-wave variety. The important question therefore is the nature of the quantum- 
critical fluctuations. Only a limited set of experiments are at present available on good 
single crystals to answer this question and no conclusive statements can be made yet. 

Two experiments suggest that quantum-criticality may be in the same universality class 
as the Cuprates, i.e. the fluctuations have a weak low frequency singularity and a broad 
nearly constant distribution in frequency. One is the measurement of thermopower for hole 
doped compound to deduce the electronic contribution to entropy in the normal state, see 



fig. (16) and the other is the resistivity measurements, also shown in the same figure but for 
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FIG. 16: On the left: Thermopower [95] in K-doped FeAs2 indicating a possible TlnT electronic 
entropy near the doping for highest T c , indicating possible quantum-criticality; On the right: Re- 
sistivity [90] in Co-doped FeAs2 for various concentration, showing a possible linear in T resistivity 
near the doping of highest T c , indicating scattering from a singular fluctuations spectra consis- 
tent with that which gives TlnT for the electronic entropy. Both are charaterisitic of Marginal 
Fermi- liquids. 

the hole-doped compound. The thermopower measurements are consistent with an entropy 
oc T ln(T) down to T c from an upper cut-off of about the room temperature in a region close 
to the highest T c . Similarly, the resistivity is linear in temperature down to T c in a similar 
region. Neither of these properties are characteristic of fluctuations around a Gaussian 
quantum critical point, especially in 3d systems [31] . 

More experiments are needed to investigate whether the Fe-Pnictides have the same 
class of quantum-criticality as the Cuprates and the heavy-fermions. Experiments which 
would be helpful are measurements of the Raman and inelastic neutron scattering spectra, 
measurements of single-particle scattering rates by ARPES measurements and of transport 
scattering rate through optical conductivity measurements. Also analysis of ARPES spectra 
in the superconducting and normal state following such experiments to decipher the spec- 
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trum of the pairing glue will undoubtedly be forthcoming. Inelastic scattering in the normal 
state to indicate the variation of the spectra with (q, u) to see whether the bulk of it has 
very slow q-dependence and a high energy cut-off as well as oj/T scaling would be especially 
helpful. 

On the theoretical issues, several models were proposed which are misguided multi-orbital 
generalizations of models for Cuprates. The opinion is converging to models similar to 
those with Hund's rule couplings and local repulsions of magnitude similar to the band- 
width. This is the class of models relevant to itinerant magnetism in metals like Ni 
or Fe. For such models, spin-polarized variational band-structure models give the right 
values for ground state magnetization and fermi-surfaces. Solution of such models for their 
fluctuation spectra is quite another matter. The physics of itinerant ant-iferromagnetism 
and ferromagnetism in such situation are not fermi-surface effects, the spin-gaps are more 
than an order of magnitude larger than the transition temperatures and essentially the entire 
band is affected. Possible large spin-gaps in the AFM pncitides should be looked for in spin- 
polarized photoemission. This always means that critical modes are soft over a large range 
in q-space leading to small intrinsic correlation lengths and huge reductions of transition 
temperatures due to the large entropy of the collective modes. There are no really good 
theories of this even for the classical fluctuations above finite temperature phase transitions. 
How this goes over to zero-temperature transitions where the soft modes may become scale- 
invariant over the entire momentum range is a marvelous unsolved question. Understanding 
of such quantum-criticality may be the key to understanding superconductivity in such 
systems. 

F. The case of Ba\_ x K x BiO?, 

I wish to single out the special case of Ba\- x K x BiOz (and BaBi\_ x Pb x Oz) because they 
have s-wave superconductivity which is most likely driven by e-e processes. This as explained 
below may be the key to superconductivity with much higher transition temperatures. 

As seen in fig. ([I]), the T c of Bai_ x K x BiO^ is an order of magnitude larger than that 
of other metals of similar low electronic density of states at the chemical potential. It was 
suggested [97] that the superconductivity could not be due primarily to electron-phonon 
interactions and some detailed calculations [HE] support this. Bi is one of about a dozen 
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elements in the periodic table which skip a valence in any of the compounds they form; 
the formal valence of Bi is always 3 + or 5 + , except for the anomaly of BaBiO^, where it 
is 4 + . In the cubic structure, this compound has 1 electron per unit-cell but the crystal 
structure is distorted producing a diamagnetic insulator with two inequivalent sites for Bi, 
one with an oxygen co-ordination favoring Bi 3+ and the other Bi 5+ . On doping with K, 
this distortion disappears in what appears to be a first-order fashion to a metallic phase 



with a high transition temperature. The phase diagram is shown in fig. (17). 




FIG. 17: T c vs. x in Ba\- x K x BiO?,\ the region for small x has incommensurate structural order 
with inequivalent sites for the nominally Bfi + and Bi 5+ . The diagram is taken from Ref.(99) 



The proper definition of an effective local interaction energy U(x + ) of an ion at a valence 
x + is the difference of the ionization and the affinity energy, i.e. 



U{x + ) = E((x + 1) + ) + E((x - 1)+) - 2E(x^ 



(45) 



Therefore, for example the positive affinity energy of spin-1/2 configuration Cu(2+) to 
spin-0 configuration Cu(l+) and the positive ionization energy to the spin-1 configuration 
Cu(3+) is equivalent to a repulsion U so that magnetic states are favored. Correspondingly, 
ions like Bi(A + ), which have a negative affinity and a negative ionization energy may be 
modelled by an attractive local interaction U. When the average valence is the skipped 
valence and the ions fluctuate quantum-mechanically between the two valences on either 
side, the fluctuations are that of a local Cooper pair. Superconductivity requires the further 
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process of phase coherence between such fluctuating pairs at different sites induced by the 
kinetic energy. 

The chemistry and physics leading to this "negative" U is not well understood. It doubt- 
less depends on the fact that the valence on either side of the skipped valence forms closed 
shells so that U calculated from the gas-phase ionization and affinity energy is much lower 
than for adjacent elements [97J, although still positive. This is likely to be reduced in a low 
energy model by polarization processes and by nearest neighbor repulsions [97] . But a quan- 
titative understanding of why in the solid state, irrespective of their chemical surroundings, 
such elements always skip valence needs further investigation. 



The phase diagram of the form in fig. (17) and some of its other properties were obtained 
in a model with negative U and repulsive nearest neighbor interactions [97]. A detailed 
investigation on this model has been carried out [100]. I wish to highlight only the experi- 
mental features which suggest that this compound has unusual properties, besides the high 
T c indicative of pairing due to electron-electron interactions. A limitation on our knowl- 
edge of this compound is that the extrapolated residual resistivity of the samples measured 
suggests a mean-free path due to impurities of only a few lattice constants. I hope better 
samples of the compound can be made and its properties in the pure limit studied more 
thoroughly. 

In the insulating state at low doping, transport and optical gaps differ by nearly an 
order of magnitude. In BaBiO^, the transport gap is 0.24 eV while the optical gap is 
about 2 eV |100] . This remarkable deviation has been well explained in the model of local 
electronic attraction |100j by the formation of Cooper pairs in the insulating state. The 
ac conductivity in the metallic state deviates from the Drude form and resembles in some 
experiments the form in the Cuprates. Tunneling measurements |101j show a conductance 
G(V) oc |V| for |V| » T. This is characteristic of a single particle scattering rate cx \cj\, the 
energy measured from the chemical potential. This is of the same form as observed through 
ARPES in Cuprates near 'optimal' doping, as well as with the peculiarity in the optical 
conductivity. What is even more interesting is that a systematic study of the slope (j(V)/| V| 
as a function of doping shows a correspondence with the measured T c . These results are 



shown in fig. (18). This suggests, unsurprisingly, that the fluctuations responsible for the 
anomalous normal state scattering are also responsible for pairing. 



There is evidence for anomalous fluctuations in Raman scattering, shown in Fig. (14). 



54 




LOG ZERO BIAS CONDUCTANCE (mhos) 

Volloge (millivolts) 

FIG. 18: Tunneling data from Sharif! et al. [lUlj and connection of the slope of conductance to Tc 

They have the same form as the critical continuum observed in the Cuprates. However, 
the \u\ dependence of the single-particle scattering rate and in transport requires that such 
a critical continuum exist over a large momentum region. There are no experiments to 
look for this nor is there any theory which leads from a negative \U\ model to such critical 
fluctuations. While more experiments in better samples are called for, a tentative case can be 
made on the basis of the Raman scattering and the tunneling data that the superconductivity 
in Bai- x K x BiOs is also due to the new universality class of fluctuations. 

The importance of this compound is that it has isotropic singlet superconductivity. From 
the considerations earlier, this is the class of electronically induced superconductivity which 
has the best chance of having the highest T c , because (1) the cut-off scale of the fluctuations 
is an electronic energy scale and (2) because the effects of inelastic scattering and self-energy 
are least deleterious for T c in the s-wave channel. But T c in Ba\^ x K x BiO^ is unimpressive 
because of its exceptionally low electronic density. 

That Bi is not unique in regard to skipping valence is shown experimentally by the 
example of Tl doped PbTe. Here at very low doping evidence |102j of fluctuations between 
the two valences, Tl l+ and T/ 3+ , i.e. of Cooper-pair fluctuation or charge Kondo effect 
[103], has been presented. This turns into a superconducting state with further doping as 
the kinetic energy increases. T c rises to about 1.5K at about 1.5% doping. Unfortunately, 
the compound cannot be doped further. 

Fabrication of compounds with valence skipping ions listed in ref. (97) is suggested so 
as to stabilize the skipped valence as the average valence, hopefully with large electronic 
density. Unfortunately, they generally seem to be unstable towards one or the other side 
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of the skipped valence; so special fabrication techniques may be required. This appears a 
likely route to high T c with the added benefit of being in compounds which are more or less 
isotropic. 



V. APPENDIX: RELATION AMONG PARAMETERS DETERMINING 
A IN E-PH INTERACTIONS 

1. "Constancy" of JV(0) < I 2 > 

Consider the Hamiltonian for the non-degenerate band of tight-binding orbitals 

H = ^ tijafa ja = ^e k a+ CT a kCT ; e k = nt cos(k.a). (46) 

i<j,cr k,(T 

A representation of the electron-phonon interaction which very efficiently includes the local 
field effects is a model in which the tight binding orbital moves rigidly with the ions. The 
electron-phonon interactions are then obtained by simply modulation of the transfer integral: 

dt-- 

%(r) wt ii ( a ) + -^-.(u i -u i )+... , (47) 

where u, is the displacment of the ion at site i from equilibrium. Moreover 

dt t a . . 

FT' 48 

or r |a| 

where tq is the characteristic radius of the tight binding orbital. This yields an elecron- 
phonon Hamiltonian 



H im = ^/(^^KA+^iKu + &J„), (49) 

where 

J k,q = i -f I £q (^ka - ^k+qa) • (50) 
| (X\T o 

Here v^ a is the velocity, h~ 1 dek/dk a . 

The parameter < I 2 > occurring in A is the average of the coupling function i£ q over 
the fermi-surface: 

< p >= /w^wv. (51) 



With (50), this is easily evaluated to give [15] 



N{0) < I 2 >^ (l/r 2 ) J2^ = K/rl (52) 

k<k F 
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where E® is the cohesive energy. This relation was used [15] to understand the near-constancy 
of N(0) < I 2 > in transition metals and their alloys. 

A modification of this calculation [TS] to include the first correction due to the non- 
orthogonality S of the nearest neighbor tight binding orbitals gives that 

N{0) < I 2 (1 ± S)E°Jr 2 , (53) 

where the (±) sign holds if the fermi-level falls in the antibonding/bonding part of the 
band. This reflects merely that the velocity in the antibonding/bonding parts of the band 
is larger /smaller approximately by (1 ± S). 

2. Relation of average electron-phonon scattering to average lattice stiffness 

We wish to derive [15] the observed relation shown in Fig.(?? between < I 2 > and the 
average lattice stiffness M < u> 2 >. The phonon frequencies u q are renormalized from 
their rigid-ion values Q q due to the creation of particle-hole pairs created by the lattice 
deformation: 

u 2 ^n 2 q + 2u q U(q,0) (54) 

where II (q, 0) is the self-energy 

U(q,0) = ^\g kMq \ 2 X (k,k + q) (55) 

k 

and x{k,k + q) is the electronic polar izability at zero frequency. 

The average rigid-ion stiffness is given by the second-derivative of the cohesive energy of 
the lattice: 

M <^>=E^o (56) 

a a 

Contributions to E(R) come from (1) The sum of the one-electron (band-structure energy) 
discussed above and (2) core-core repulsion energy minus the electron-electron repulsion 
energy, the latter having been counted twice in a self-consistent band-structure energy. 
For any long-wavelength property (2) sums to zero and it is negligible for consideration of 
the average stiffness compared to that of (1) because it is slowly varying with inter-site 
separation while (1) varies on the scale of the variation of the overlap integrals. In fact 
chemists calcualting the structure and vibration frequency of metals have found that good 
results are obtained by neglecting it altogether. For transition metals and compounds, we 
may note further the justification that about 60 — 80% of the cohesive energy is provided 
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by d-electrons alone. Even their crystal structure is correctly predicted on the basis of (1) 
alone. 

Taking the second derivative of E(R) 

M < n 2 >^ (1 ± S)E°Jrl w JV(0) < I 2 > . (57) 



Consider the second term in (54). If Uk,k+q is the residual electron-electron interaction 
between two tight-binding orbitals, 

X (k, k + q) = X °(k.k + q)/(l + U kMqX °(k, k + q)) (58) 

where x°{k, k + q) is the zero order electronic polarizability at zero frequency, 

rW/2 f W/2 f( N _ f ( A 

X°(k,k') = -(1/7T 2 ) / del de' J -^ J -^ImG{k,e)ImG{k',e') (59) 

J-w/2 J-w/2 e — e' 

Using (l/ir)ImG(k, e) = 5(e — ek) and ignoring logarithmic corrections, this for typical k, k'is 

just iV(0) 2 W. Therefore using Eq. (49, 50 ) for g(k,k') gives 

2^Afw q |^ )fc+ff |V(k,k + q) w -iV 2 (0)lU < I 2 > (60) 

k,q 

For most interesting cases A^(0)iy > 1, so that we must consider the change of phonon 
frequencies due to the electronic polarizability and local field corrections. For large values 
of k, k + q, which alone are of interest, we may take the effective interaction Uk y k+q ~ W, as 



one finds in the t-matrix approximation. This estimate in Eqs.(52, 54) gives 

N(Q) < I 2 > /M <u 2 > ps (1 + N(0)W) (61) 

This derivation is necessarily crude. It works qualitatively (and semi-quantitiatively with an 
overall correction factor of about 1/4) because we are calculating average of a quantity over 
the entire zone. These results are discussed in relation to Fig. ^ earlier. The basic physical 
points used, that local field corrections or modulation of bonds due to lattice vibrations are 
essential in understanding the details of e-ph interactions and phonon dispersion and that 
the parameters determining T c are inter-related, is well borne out in detailed calculations of 
phonon frequencies and electron-phonon couplings and their comparison with experiments 
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